Diversity of critical behavior within a universality class 
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We study spatial anisotropy effects on the bulk and finite-size critical behavior of the O(n) sym- 
metric anisotropic ip 4 lattice model with periodic boundary conditions in a d-dimensional hypercubic 
geometry above, at and below T c . The absence of two-scale factor universality is discussed for the 
bulk order-parameter correlation function, the bulk scattering intensity, and for several universal 
bulk amplitude relations. The anisotropy parameters are observable by scattering experiments at 
T c . For the confined system, renormalization-group theory within the minimal subtraction scheme 
at fixed dimension d for 2 < d < 4 is employed. In contrast to the e = 4 — d expansion, the fixed-d 
finite-size approach keeps the exponential form of the order-parameter distribution function unex- 
panded. For the case of cubic symmetry and for n = 1 our perturbation approach yields excellent 
agreement with the Monte Carlo (MC) data for the finite-size amplitude of the free energy of the 
three-dimensional Ising model at T c by Mon [Phys. Rev. Lett. 54, 2671 (1985)]. The e expansion 
result is in less good agreement. Below T c a minimum of the scaling function of the excess free en- 
ergy is found. We predict a measurable dependence of this minimum on the anisotropy parameters. 
The relative anisotropy effect on the free energy is predicted to be significantly larger than that on 
the Binder cumulant. Our theory agrees quantitatively with the non-monotonic dependence of the 
Binder cumulant on the ferromagnetic next-nearest neighbor (NNN) coupling of the two-dimensional 
Ising model found by MC simulations of Selke and Shchur [J. Phys. A 38, L739 (2005)]. Our theory 
also predicts a non-monotonic dependence for small values of the antiferromagnetic NNN coupling 
and the existence of a Lifshitz point at a larger value of this coupling. The nonuniversal anisotropy 
effects in the finite-size scaling regime are predicted to satisfy a kind of restricted universality. The 
tails of the large-L behavior at T / T c violate both finite-size scaling and universality even for 
isotropic systems as they depend on the bare four-point coupling of the <p 4 theory, on the cutoff 
procedure, and on subleading long-range interactions. 



I. INTRODUCTION AND SUMMARY 

A major achievement of the renormalization-group 
(RG) theory is the proof that critical phenomena can 
be divided into distinct universality classes (for a review 
see, e.g., [l|). They are characterized by the spatial di- 
mension d and the symmetry of the ordered state which, 
for simplicity, we assume in the following to be O(n) sym- 
metric with an n component order parameter. (For other 
universality classes see, e.g., Q.) Within a given (d,n) 
universality class, all bulk systems (with finite-range in- 
teractions and with subleading long-range interactions of 
the van der Waals type) have the same critical exponents 
and the same thermodynamic functions near criticality in 
terms of universal scaling functions that are obtained af- 
ter a rescaling of two amplitudes : that of the singular 
part of the bulk free energy density f Sy b and that of the 
field h conjugate to the order parameter. This is sum- 
marized in the asymptotic (small t — (T — T c )/T c , small 
h) scaling form (below d — 4 dimensions) 



f a , b (t, h) = A^tr W ± (A 2 h\t 



-ps\ 



(1.1) 



with universal critical exponents v, [3, 5 and the univer- 
sal scaling function W±(z) above (+) and below (— ) T c . 
Once the universal quantities are known one knows the 
asymptotic thermodynamic critical behavior of all mem- 
bers of the universality class provided that only the two 
nonuniversal amplitudes A\ and A 2 are specified. (For 
the application to real systems, additional experimental 
information is necessary to identify the order parame- 



ter and the appropriate thermodynamic path tangential 
to the coexistence line.) We refer to this property as 
thermodynamic two-scale factor universality. Here uni- 
versality means the independence of all microscopic de- 
tails such as lattice structure, lattice spacing, and the 
specific form and magnitude of the finite-range or sub- 
leading long-range interaction. This implies that both 
fluids and anisotropic solids within the same universality 
class have the same scaling function W±. 

This important concept of scaling and thermodynamic 
two-scale factor universality was extended to the distance 
(r) dependence of bulk correlation functions Q and to 
the size (L) dependence of quantities of confined systems 
[HIE El (for reviews see, e.g., @, HI)- It is this extended 
hypothesis which is in the focus of the present paper. We 
shall present results for the finite-size critical behavior of 
the free energy above, at and below T c that demonstrate 
a considerable degree of diversity within a given (rf, n) 
universality class primarily due to spatial anisotropy in 
lattice systems with non-cubic symmetry, but also due 
to the lattice spacing a in systems with cubic symme- 
try and due to the bare four-point coupling u of the 
(p theory even in the isotropic case. In this context we 
also discuss nonuniversal effects related to the cutoff and 
to subleading long-range (van der Waals type) interac- 
tions. This diversity suggests to distinguish subclasses of 
interactions within a given universality class where the 
subclasses have different bulk amplitude relations, differ- 
ent bulk correlation functions, and, for given geometry 
and boundary conditions (b.c), different finite-size seal- 
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TABLE I: Subclasses of asymptotic critical behavior within a (d, n) universality class for O(n) symmetric systems in a cubic 
volume V — L d with periodic boundary conditions for general n above T c and n — 1 below T c . All subclasses have the same fixed 
point value u* (d, n) of the renormalized four-point coupling, the same critical exponents, and the same bulk thermodynamic 
scaling functions. This table complements Table IV of [J]. 



classes 
of interactions 
6K(k) in P7) 


basic lengths, 
nonuniversal parameters 


bulk 
amplitude 
relations 


bulk 

correlation functions 


finite-size 
effects 


isotropic short 
range" 
k 2 + 0(k 4 ) 


correlation length 6 £±, 
two nonuniversal 
amplitudes Ci , C2 , 
four-point coupling Mo 


two-scale 
factor 
universality 


r /£± ^ O(l) : universal 
isotropic power-law scaling form; 
r £± : exponential form with 
nonuniversal tails 


L/f± < O(l) : universal 
power- law scaling form; 
L 2> f± : exponential form 
with nonuniversal tails 


anisotropic short 
range c 

det A > 

A a @ = A( ja 


d principal correlation 

lengths f , 
up to l)/2 + 1 
nonuniversal parameters'* 
Cii C 2 , A a) 3, 
four-point coupling Mq 


multi- 
parameter 
universality 


r /£± ~ O(l) : universal 
power-law scaling form with 
d(d + l)/2 + 1 nonuniversal para- 
meters in the scaling arguments; 
r 2> £± : exponential form with 
nonuniversal tails 


L/£jf^ < O(l) : nonuniversal 
power- law scaling lorm, 
restricted universality; 
L S> f± : exponential form 
with nonuniversal tails 


isotropic 
subleading long 
range 6 
k 2 - b |k| CT 
2 < a < 4 


correlation length £±, 
interaction length scale 

6 V(-2) ) 

five nonuniversal 
parameters Ci , C2, 6, <x, uq 


two-scale 
factor 
universality 


r /£± i; O(l) : universal 
power-law scaling form; 
r/£± > O(l) : 
nonuniversal power-law form 
depending on &, a 


L/£± < O(l) : universal 
power-law scaling form; 
L/C±>0(1): 
nonuniversal power-law form 
depending on b, a 



° Refs. hie n, [is [ni 

6 For isotropic systems, and £_ denote the second-moment 
bulk correlation lengths above and below T c , with a universal ra- 
tio For anisotropic systems, £v"' are the principal bulk 

correlation l eng ths with universal ratios fj ''/?— > a = ^> % •••> 
c Refs. [H, [T3l 

d The reduced anisotropy matrix A = A/(det A) 1//d has d(d + 
l)/2 — 1 independent matrix elements A a g. 

e Refs. [HI, [H, m, EE E3, Ell 



ing functions. All of these nonuniversal differences exist 
in the asymptotic critical region \t\ <C 1, L 3> a, and 
r>o where corrections to scaling in the sense of Weg- 
ner fl9j | are negligible. A summary of these properties is 
given in Table I. The basic framework of RG theory is 
fully compatible with this diversity of critical behavior. 

Spatially anisotropic systems such as magnetic mate- 
rials, alloys, superconductors |20l|. and solids with struc- 
tural phase transitions [2TL |22| represent an important 
class of systems with cooperative phenomena. One may 
distinguish between long-range anisotropic interactions 
(such as dipolar, RKKY, and effective elastic interac- 
tions) and short-range anisotropic interactions which in- 
clude the Dzyaloshinskii-Moriya-type antisymmetric ex- 
change [23[ and the spatially anisotropic Heisenberg ex- 
change interactions which, in the long- wavelength limit, 
are described by a d x d anisotropy matrix A [I2I [l3| . 
We shall confine ourselves to a detailed study of the lat- 
ter type of systems but the general aspects of our results 
have an impact also on the former type of anisotropic 
systems and on systems of other universality classes [2j, 
for example on the range of validity of universality for 
anisotropic spin glas ses [24| or for anisotropic surface crit- 
ical phenomena [2a ]. 

A characteristic feature of spatial anisotropy with non- 



cubic symmetry is the fact that there exists no unique 
bulk correlation-length amplitude but rather d different 
amplitudes £q™' in the directions a = 1, d of the d prin- 
cipal axes. Such systems still have a single correlation- 
length exponent v provided that det A > 0. (We do not 
consider strongly anisotropic systems with critical expo- 
nents different from those of the usual (d, n) universality 
classes, see e.g. [26|].) Non-cubic anisotropy effects in 
crystals with cubic symmetry can be easily generated by 
applying shear forces. In earlier work on two-scale fac- 
tor universality [H, H, H3, [H, [29|, [3(|, isotropic systems 
with a single bulk correlation length ^ were considered 
and important universal bulk amplitude relations were 
derived that depend on only two nonuniversal parame- 
ters. Recently some of these relations were reformulated 
for anisotropic systems within the same universality class 
[T2I [l3j . In Sect. Ill of the present paper we give a 
derivation of these and other relations above and below 
T c and express them in terms of universal scaling func- 
tions. The physical quantities entering these relations 
depend, in general, on d{d + l)/2 + 1 nonuniversal pa- 
rameters. We also present the appropriate formulation of 
the bulk scattering intensity of anisotropic systems near 
criticality in terms of the eigenvalues of the anisotropy 
matrix and discuss the nonuniversal properties of bulk 
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correlation functions. 

For confined systems with a characteristic length L the 
hypothesis of two-scale factor universality is summarized 
by the asymptotic (large L, small t, small h) scaling form 
for the singular part of the free energy density (divided 
by k B T) 0,0,1 

f 8 (t, h, L) = L- d TidtL 1 '", G 2 hL^l v ). (1.2) 

where T(x,y), for given geometry and b.c, is a univer- 
sal scaling function and where the two constants C\ and 
C 2 are universally related to the bulk constants A\ and 
A2 of p.l[) . For simplicity we shall confine ourselves to 
a hypercubic shape with volume V = L d and with pe- 
riodic b.c. Calculations of f s (t,0,L) for this case were 
carried out within the spherical model [3l[ which sup- 
ported the form of (|1.2p . For n = 1 the scaling form 
(|1.2|) was discussed in the framework of the e = 4 — d ex- 
pansion [33, [33[ ■ No theoretical prediction for the func- 
tion T(x, y) is available up to now for finite n in cubic 
geometry, except in the large-n limit (TlJ. Monte Carlo 
(MC) simulations [H, Hi| [36| for three-dimensional Ising 
models with nearest-neighbor (NN) couplings on different 
lattices of cubic symmetry were consistent with the uni- 
versality of the amplitude ^(0,0). These models belong 
to the subclass of (asymptotically) isotropic systems. 

It was already noted in 0,13,1221 that lattice anisotropy 
is a marginal perturbation in the RG sense, thus it was 
not obvious a priori to what extent two-scale factor uni- 
versality is valid in the presence of anisotropic couplings 
Q. It was also known that, for most anisotropic systems, 
(asymptotic) isotropy can be restored by an anisotropic 
scale transformation [H, [3!| (for further references see 
[l2l]). Recently it was pointed out [12] that, in systems 
with anisotropic interactions of non-cubic symmetry, the 
scaling function T is indeed affected by anisotropy. In 
particular, it was shown [l3| that by means of an ap- 
propriate rescaling of lengths a transformation to an 
(asymptotically) isotropic system is always possible pro- 
vided that the anisotropy matrix A is positive definite 
and that the rescaling is performed along the d nonuni- 
versal directions of the principal axes which, in general, 
differ from the symmetry axes of the system. This rescal- 
ing is equivalent to a shear transformation which distorts 
the shape, the lattice structure, and the boundary condi- 
tions in a nonuniversal way (e.g. from a cube to a paral- 
lelepiped, from an orthorhombic to a triclinic lattice, and 
from periodic b.c. in rectangular directions to those in 
non-rectangular directions). This nonuniversality is re- 
flected in a dependence of the scaling function T on the 
anisotropy matrix A, in addition to the dependence on 
C\ and C2. 

Specifically, on the basis of the results of renormalized 
perturbation theory in Sects. IV - VI, we propose that, 
for anisotropic systems with the shape of a cube, (|1.2p is 
to be replaced by (4p| 



with V = L(det A)- 1 /^), ti =_h(det A) 1 / 4 , and 
with the reduced anisotropy matrix A = A/(det A) 1 /^, 
det A > 0. The nonuniversal constants C[ and C' 2 will 
be specified in Sect. VI in terms of the asymptotic am- 
plitudes £q + and £' c of the second-moment bulk corre- 
lation lengths for T > T c , h' = and for T = T c , 
h! 7^ 0, respectively, of the transformed isotropic sys- 
tem. The free energy density f' s = /^detA)" 1 / 2 of the 
parallelepiped with the volume V = V(det A) -1 / 2 and 
with A' = A'/(det A') = 1 (isotropy) then attains the 
scaling form 

f' s (t,h',L') = L'- d ^ A (C[tL fl /",C' 2 h'L'^) (1.4) 

where the characteristic length L' — V n ^ d determines 
the overall size of the parallelepiped and 



^ iso .a{x,v) = Fcube(x,y, A) . 



(1.5) 



Equation (|1.4[) has the structure of the isotropic 
Privman-Fisher scaling form (|1.2p with a rescaled length 
L' and with only two nonuniversal constants C[ and C' 2 
which, superficially, appears to be in agreement with two- 
scale factor universality. The remaining d(d + l)/2 — 1 
nonuniversal parameters, however, are hidden in the in- 
dex "iso, A". This index is the notation for a sys- 
tem with the shape of a parallelepiped whose interaction 
5K'(\s!) = k' 2 + 0{k' A ) is (asymptotically) isotropic and 
whose d(d — l)/2 angles and d— 1 length ratios are deter- 
mined by the d(d + l)/2 — 1 nonuniversal parameters of 
the reduced anisotropy matrix A. These parameters ap- 
pear in the calculation of T iso ^ via the summation over 
the discrete k' vectors in the Fourier space of the paral- 
lelepiped system since the k' vectors depend explicitly on 
A, unlike the k vectors of the cubic system. Thus for the 
calculation of T iso ^ the same nonuniversal information 
is required as for the calculation of T cu b e {x, y; A). 

For general A the function T cu b e {x, 0; A) was pre- 
sented in |l2( for t > in the large-n limit. Furthermore, 
quantitative predictions were made for the nonuniversal 
dependence of the critical Binder cumulant d, [cO ] 



U C ube (-^-) 



d 4 T cube (0,y;A)/dy 4 



(9 2 ^ cube (0,y;A)/5y 2 



(1.6) 



for n = 1,2,3 both in threejlljll] and two [H di 



f,(t,h,L) = L- d T cube {C[tL lx l v ,C' 2 h'L'^' v -K) 



(1.3) 



mensions. MC simulations [42|, |43j for the anisotropic 
three-dimensional Ising model indeed showed nonuniver- 
sal anisotropy effects which, however, did not agree with 
the theoretical prediction. More accurate MC simula- 
tions [44| for the anisotropic two-dimensional Ising model 
demonstrated the nonuniversality of the critical Binder 
cumulant but no comparison with a quantitative theo- 
retical prediction was available for this two-dimensional 
case. Thus the anisotropic finite-size theory of [HI, [l3| is 
as yet unconfirmed. 

In Sects. IV - VI of this paper we derive the finite- 
size scaling function F^ be {x, 0; A) of the singular part of 
the excess free energy density ff x — f s — f s>b at h = 



4 



above, at, and below T c for the n = 1 universality class 
in 2 < d < 4 dimensions on the basis of the anisotropic 
ip 4 lattice model. For the isotropic case at T c we find 
excellent agreement with the MC data of Mon [U H|| . 
Slightly below T c we find a minimum of the scaling func- 
tion that is similar to the minimum of the scaling func- 
tion of the critical Casimir force for the d = 3 Ising 
model in slab geometry with periodic boundary condi- 
tions (H, [46j. For future tests of our theory by MC 
simulations we consider both three- and two-dimensional 
anisotropics. In both cases we predict a measurable de- 
pendence of the minimum on the anisotropy parameters, 
thus demonstrating the nonuniversality of the finite-size 
scaling function of the excess free energy density. The 
magnitude of this anisotropy effect is predicted to be con- 
siderably larger than that on the Binder cumulant. 

We believe that a similar nonuniversal dependence can 
be derived for the critical Casimir force by means of our 
perturbation approach. For n — ► oo and at T = T c , 
the nonuniversality of the Casimir amplitude due to 
anisotropy was already demonstrated within the <p 4 the- 
ory in [12j. This suggests that, for given geometry and b. 
c, the existing theoretical 0, B El, Hfl, HH \M, H HI 
and MC [H, [5| [H, HH results for the Casimir force scal- 
ing function are not universal within the entire univer- 
sality class but are restricted to the subclass of isotropic 
systems. Extensions to the subclass of anisotropic sys- 
tems are, in general, not straightforward and cannot 
be obtained just by transformations but require new 
nonuniversal input, new analytical and numerical calcu- 
lations, and new MC simulations. Experiments, e.g., in 
anisotropic superconducting films pol . [56j , could, in prin- 
ciple, demonstrate the nonuniversality in real systems. 

Our present results for the p 4 theory cannot be applied 
directly to two-dimensional critical phenomena. Never- 
theless we are able to study two-dimensional anisotropy 
effects within a three-dimensional model. For the pur- 
pose of a comparison with the two-dimensional MC data 
[44j |. we consider (in Sect. VIII) a three-dimensional <p 4 
lattice model with the same two-dimensional anisotropy 
in the horizontal planes as in the two-dimensional model 
Ising model studied by Selke and Shchur [44| . Our theory 
agrees quantitatively with the non-monotonicity of the 
Binder cumulant as a function of the anisotropiy ferro- 
magnetic next-nearest neighbor (NNN) coupling found in 
[ii]]. We also predict a non-monotonicity for small anti- 
ferromagnetic couplings and the existence of a Lifshitz 
point at a larger value of this coupling. Very recent pre- 
liminary MC data by Selke [|37| for the too-dimensional 
Ising model indeed reveal such a non-monotonicity that 
was not yet detected in (44J. We predict a similar 
anisotropy effect for the excess free energy density of 
the anisotropic two-dimensional Ising model. This ef- 
fect can become quite large if one of the eigenvalues of 
A approaches zero, in particular if a Lifshitz point is ap- 
proached (Sect. VIII). 

An important property of the scaling form (|1.3|) is that 
it depends on A but not on other nonuniversal param- 



eters such as the bare four-point coupling, the lattice 
spacing, and the cutoff of ip 4 field theory. This is a kind 
of restricted universality since it implies that the same 
finite-size scaling functions exist for the large variety of 
those systems within a universality class that have the 
same reduced anisotropy matrix A (and the same ge- 
ometry and boundary conditions). In Sect. IX we pro- 
pose two examples for testing this hypothesis of restricted 
finite-size universality by MC simulations for spin models 
with anisotropic interactions. For recent tests of finite- 
size universality of two-dimensional Ising models with 
(asymptotically) isotropic interactions see [58|, [5^| (see 
also Table 10.1 of 0). 

Unlike the bulk scaling function W±(z), (II. ip . that 
is valid in the entire range — 00 < z < 00 of the scal- 
ing argument z, the finite-size scaling functions such as 
^cube(x, y \ A) are valid only in a limited range of x and y, 
above the shaded region in Fig. 1. In the shaded region, 
nonuniversal nonscaling effects become nonnegligible and 
even dominant for sufficiently large |x| and \y\ for both 
short-range and subleading long-range interactions. In 
this region, not only the correlation lengths are relevant 
but also additional nonuniversal length scales such as the 
lattice spacing a, the inverse cutoff A -1 , the length scale 
Uq 1 ^ set by the four-point coupling , and the van-der- 
Waals interaction-length oVO- 2 )^ as discussed in Sect. 
X. 

A (l/L') 1 /" 



lowest - modt separation , 
N N perturbation approach ' 

ordinary ^ / ordinary 



perturbation N | ' perturbation 



approach N \ 


nonuniversal 


/ approach 




Kb 






scaling 




nonscaling 




nonscaling 







t 



FIG. 1: Asymptotic part of the L'^ 1 ^ — t plane at ft = for 
the anisotropic if 4 theory in a cubic geometry with periodic 
boundary conditions. In the central finite-size region (above 
the dashed lines), the lowest mode must be separated whereas 
outside this region ordinary perturbation theory is applicable. 
Above the shaded region, finite-size scaling is valid but with 
scaling functions that depend on the anisotropy parameters 
A a p, see (|1.3fl . In the large - L' regime at t 7^ (shaded re- 
gion) finite-size scaling and universality are violated for both 
short-range and subleading long-range interactions and for 
both isotropic and anisotropic systems. A similar plot is valid 
for the L'- 0S/ " - ti plane at T = T c . 

We briefly comment on the methodological aspects of 
our finite-size calculations. As far as the field-theoretic 
[60| renormalization of the ip lattice model is concerned 
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we use the minimal subtraction scheme [6l| not within 
the e expansion but at fixed dimension d, as introduced 
in [HH and further developed in [63[. As far as finite-size 
theory is concerned we further develop earlier approaches 
[32l loll . [6f| [66| that have been successfully used to solve 
several finite-size problems in the past [6?], El, [§3, l7Cj |. 
After the transformation from the anisotropic to an 
isotropic system, the same renormalization constants {Z 
factors) and the same fixed-point value u* of the renor- 
malized four-point coupling are obtained as for the stan- 
dard isotropic tp 4 Hamiltonian. For this reason, the same 
fixed-point Hamiltonian and the same critical exponents 
govern isotropic and (weakly) anisotropic systems - they 
belong to the same universality class. The crucial point, 
however, is that not only the fixed-point value u* but 
also the orientation of the eigenvectors (principal axes) 
of the fixed-point Hamiltonian relative to the orientation 
of the given boundaries of the confined anisotropic sys- 
tem determine the finite-size scaling functions. This is a 
physical fact that introduces a source of nonunivcrsality 
that cannot be eliminated by transformations and that 
makes anisotropic confined systems distinctly different 
from isotropic confined systems within the same univer- 
sality class. 

The main result for f% x will be obtained in the central 
finite-size regime (above the dashed lines in Fig. 1) where 
the finite-size effects are most significant and where it is 
necessary to separate the lowest-mode from the higher 
modes. In this regime finite-size scaling is valid in the 
form of (|1.3p . We compare the result of our fixed - d 
perturbation approach [6^, [H, [H, [6(| with that of the 
e expansion approach. The advantage of the former ap- 
proach is that it keeps the exponential structure of the 
order-parameter distribution function unexpanded. This 
leads to a result at T c in excellent agreement with the 
MC data in the isotropic case [3J, |3f| and lends credi- 
bility also to the quantitative features of our predictions 
of anisotropy effects. The e expansion result at T c turns 
out to be in less good agreement. 

The separation of the lowest mode is inadequate in 
the limit of large L' 3> £'± at fixed T ^ T c . In order to 
capture the exponential structure of finite-size effects for 
large V we complement (in Sect. X) our results by or- 
dinary perturbation theory outside the central finite-size 
regime (below the dashed lines in Fig. 1). This includes 
a small but finite region where finite-size scaling is vio- 
lated (shaded region in Fig. 1). There exists diversity 
rather than universality of finite-size critical behavior in 
this region depending on all microscopic details of the in- 
teractions such as the lattice spacing, the bare four-point 
coupling, the cutoff of the ip 4 theory, and the amplitude 
of subleading long-range interactions. This diversity can 
be traced back to a corresponding diversity of the large- 
distance (r' 3> £' ± ) behavior of the bulk order parameter 
correlation function Gt, [Io| where r' is the distance in 
the transformed isotropic bulk system, as discussed in 
Sect. III. For Gb there exists a region of the r' _1 — ^^T 1 
plane (shaded region in Fig. 2) that is the analogue of 



n 1/r' 



nonuniversal 

scaling 
arguments 



FIG. 2: Asymptotic part of the r' _1 — ^ 1 plane at h = for 
anisotropic bulk systems. Above the shaded region, there ex- 
ists a universal scaling function $±(r' /£'±,0) of the bulk cor- 
relation function Gt, (|3.19|l . The scaling argument, however, 
contains the spatial variable r' = |x'|, (|3.23[1 . that depends 
on the anisotropy matrix (A a p) with d(d+ l)/2 nonuniversal 
parameters. In the large - r' regime at t 7^ (shaded region), 
scaling and universality are violated for both short-range and 
subleading long-range interactions and for both isotropic and 
anisotropic systems. A similar plot is valid in the r' _1 — 
plane at T = T c ,ti / 0, with & = i' c \h'\- u/{riS) . 



the shaded region of Fig. 1 . In the isotropic case, this re- 
gion is of physical relevance for fluids with van der Waals 
interactions [TO, 0, E EE GJ, S3, HI • 



II. ANISOTROPIC ip* LATTICE MODEL 
A. Hamiltonian with spatial anisotropy 

We start from the O(n) symmetric <p lattice Hamilto- 
nian (divided by ksT) 



H = v 



N 
i=l 



N 



(2.1) 



r (T) = r 0c + a t, t = (T - T c )/T c with a > 0, u > 0. 
The variables (pi = </?(xj) are n-component vectors on N 
lattice points Xj = (2^1, £i2, ■ ■ ■ , ^id) of a d-dimensional 
Bravais lattice with the finite volume V — Nv with the 
characteristic length L — V x l d where v is the volume of 

the primitive cell. The components ip\ , fi — 1,2, ... ,n 
of tpi vary in the continuous range — 00 < ip^ < 00. The 
couplings Ki.j — Kj^ = K(xi — Xj) and the temperature 
variable ro(T) have the dimension of L~ 2 whereas the 
variables ipi have the dimension of i( 2 ~ d )/ 2 such that H 
is dimensionless. The free energy per unit volume divided 
by k B T is 



f(t,h,L) = -V~ x hxZ 



(2.2) 
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Z(t,h,L) 



N 

n 



,n(2-d)/(2d) 



cxp(-H) (2.3) 



where Z is the dimensionless partition function. The 
total excess free energy density is defined as 



f ex (t,h,L) = f(t,h,L)-f b (t, h) 



(2.4) 



where fb(t, h) = lim^c^ /(i, h, L) is the bulk free energy 
density. Following [1,0, Q we shall decompose /, for large 
L, into singular and non-singular parts 



f(t,h,L) = f s (t,h,L)+f ns (t,L) 



(2.5) 



where f n s(t, L) has a regular t dependence around t = 0. 
In earlier work on finite-size effects it was supposed 
[1,0, EH that, for periodic boundary conditions, one can 
assume that there exists no L dependence of the non- 
singular part f ns . Adopting this assumption leads to a 
misinterpretation [ll| of the singular part f s of the free 
energy density and of the Casimir force in the presence 
of a sharp cutoff of tp 4 field theory. Here we shall not 
exclude the possibility of an L dependent nonsingular 
part f ns (t, L) even for periodic boundary conditions if 
long-range correlations are present . As will be shown in 
Sect. X, this will reconcile the earlier results [ll| with 
the concepts of finite-size scaling. 

For periodic b.c, the Fourier representations are 
ip(xj) = V- 1 Y^k e ik ' x ^(k) and 

K(xt - xy) = N- 1 e lk ' (Xl_Xj) A(k) , (2.6) 

k 

where the summations run over the N discrete vec- 
tors k of the first Brioulluin zone of the reciprocal lat- 
tice. We assume finite-range interactions Kij with a 
finite value K(0) = . =1 Kij. In terms of the 

Fourier components the Hamiltonian reads 

H = V- 1 Y, \[tq + <f£(k)]0(k)0(-k) - kip(0) 

k 

+ u V- 3 £[0(k)0(p)][0(q)0(-k - p - q)] (2.7) 

kpq 

where 5K(k) = 2[K(Q) — lf(k)]. In perturbation theory, 
ro + 5K(k) plays the role of an inverse propagator. 

The Hamiltonian H is isotropic in the vector space 
of the n-component variables ipi and <y3(k) but may be 
anisotropic in real space and k space. A variety of 
anisotropies may arise both through the lattice struc- 
ture and through the couplings Kij. They manifest 
themselves on macroscopic length scales via the d x d 
anisotropy matrix A = (A a p) and the anisotropy tensor 
B = (B a p^g) that appear in the long- wavelength form 

d d 
ex, 0=1 a,/?, 7, 5 

+ <3(fc 6 ) . (2.8) 



Odd powers of k a are excluded because of inversion sym- 
metry of the Bravais lattice. (For the case of non-Bravais 
lattices see the discussion in section II. C of Ref. (zH .) For 
cubic symmetry, A has the isotropic form A a p — cq d a fj 
while the (9(fc 4 ) terms of cubic systems differ from those 
of isotropic systems. In Sects. Ill and X we shall con- 
sider the model (|2.7j) also in a fully isotropic form with 
the short-range interaction 6K (k) = k 2 including a finite 
cutoff A and, for n = 1, with the subleading long-range 
interaction [ll|, [3, QJ, M, E3, EE M, H 



5K(k) 



Mkl' 



0(k 4 ) , 2 < a < 4 (2.9) 



with b > 0. The second term of the interaction (|2.9|) is 
usually classified as "irrelevant" in the renormalization- 
group sense [HO] since it leaves some (but not all) of the 
universal quantities unchanged: critical exponents and 
bulk thermodynamic scaling functions. This terminology 
is somewhat misleading as the term —b |k| CT changes not 
only the leading finite-size critical behavior at T ^ T c (in 
the shaded region of Fig.l) but it also destroys the uni- 
versality of the leading bulk critical behavior of the order- 
parameter correlation function Gb (and of other bulk cor- 
relation functions): Gf, attains an interaction-dependent 
power-law structure [3, Ell in the large- distance regime 
at T 7^ T c (in the shaded region of Fig. 2) whereas sys- 
tems with purely short-range interaction in the same uni- 
versality class have an exponentially decaying Gb in this 
regime (this decay has, in addition, a nonuniversal expo- 
nential tail, see Sect. X). 

The expression for A a p and B a p~ t $ in terms of the mi- 
croscopic couplings Ki j is given by the second moments 

El 



A a B = Ap a = N- 



N 



i,j=l 



(2.10) 

and the fourth-order moments of AT^j , respectively. They 
have been classified and studied in the context of the bulk 
correlation function in Ref. [Til ] . The symmetric ma- 
trix A depends only on the lattice structure and on the 
pair interactions Kij and is independent of the boundary 
conditions and the geometry of the system. Its eigenval- 
ues A Q , a = 1,2, ...,d, and eigenvectors e^"^ are deter- 
mined by the eigenvalue equation Ae'"' = A Q e(") with 
e ( Q ) . e (/3) — § a/3 ^ j n orc ier to have an ordinary critical 
point of the usual (d, n) universality classes we assume 
positive eigenvalues A Q ,det A = n i=1 A Q > 0, and that 
the fourth-order moments B a fj^s enter only the correc- 
tions to scaling. The critical point occurs at h = and 
at T — T c corresponding to some critical value ro(T c ) = 
roc that is defined implicitly by limt^o+ Xb(t, 0) = 
where Xb(t>h) — —HmL~>ood 2 f(t,h,L)/dh 2 is the bulk 
susceptibility for t > 0. This implies that r (T c ) = 
roc{Kij,v,Ua) depends on the lattice structure, on v, on 
Mo, and on all couplings Ki j. 

The matrix A affects the observable bulk critical be- 
havior : the eigenvalues X a enter the amplitudes of the 
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bulk correlation lengths in the direction of the prin- 
cipal axes; the latter are determined by the eigenvectors 
e( Q ' of A which provide the reference axes for the spa- 
tial dependence of the anisotropic bulk order-parameter 
correlation function 

Gb(xi -x.y,t,h) = lim {< tpiipj > - < ip > 2 } (2.11) 

V— *oo 

where < p > (t,h,L) = —df(t,h,L)/dh. Correspond- 
ingly, the matrix A determines the anisotropy of the k 
dependence of the Fourier transform 

G b (k; t,h)=vJ2 e- lk x G b (x; t, h) (2.12) 



which is proportional to the observable scattering inten- 
sity. The principal axes must be distinguished from the 
symmetry axes of the Bravais lattice. The latter depend 
only on the lattice points Xj but not on the couplings 
Kij . Below an example is given where the principal axes 
differ from the symmetry axes. 

The long- wavelength approximation takes into account 
only the leading 0(k a kp) term of 5K(k). In real space 
this is equivalent to using the p A continuum Hamiltonian 
for the vector field <p(x) 



H 



field 



E 



A ap dp dp 



% 2 , 

2 ^ 2 dx a dxp 

a, 0—1 K 



+ u (p 2 ) 2 - hp 



(2.13) 



with some cutoff A. 

Various types of anisotropies may result not only from 
pair interactions on rectangular lattice structures but 
also from nonrectangular lattice structures, from effec- 
tive many-body interactions as well as from distortions 
of the lattice structure, e.g., due to external shear forces. 
The semi-macroscopic continuum model (|2.13[) is ex- 
pected to be of general significance in that it provides 
a complete long-wavelength description of a large class 
of real systems near criticality whose nonuniversal prop- 
erties can be condensed into the d(d + 1)/2 parameters of 
the anisotropy matrix A, in addition to the nonuniversal 
parameters ro,UQ, h, A. The quantities A a p depend on 
all microscopic details (lattice structure, electronic struc- 
ture, many-body interactions) which, in general, are not 
known a priori for a given material. Thus the matrix 
elements A a p represent phenomenological parameters of 
a truly nonuniversal character. Consequently, physical 
quantities depending on A a p ( such as F C ube(0, 0; A), the 
Binder cumulant U cu be(A-), and the critical Casimir am- 
plitude) are nonuniversal as well. 

For an appropriate formulation of the bulk order- 
parameter correlation function (see Sect. Ill) and of 
finite-size scaling functions (see Sect. VI) it will be im- 
portant to employ the reduced anisotropy matrix A = 
A/(detA) 1 / d , Ae( Q ) = with the eigenvalues 



A Q = A a /(det Af/ d > and with det A = rj a =i X a = 1. 
The matrix A is independent of the kind of variables pi 
on the lattice points and independent of the number n of 
components of pi. It is well defined, e.g., also for mod- 
els with fixed-length spin variables Si with |Si| = 1 and 
for Ising models with discrete spin variables o~i = ±1 in- 
stead of the continuous vector variables pi. Thus the 
XY and Ising Hamiltonians Hxy = — j Ji.jSi ' Sj 
and Hi S i ng = —J2ijJi,j (J i' J j have the same reduced 
anisotropy matrix A and the same reduced eigenvalues 
A a as the p A lattice Hamiltonian if these models are de- 
fined on the same lattice points Xj and if the couplings 
Ji.j are proportional to K+j. 

As an illustration we consider an L x L x L simple- 
cubic lattice model (Fig. 3) with a lattice constant a 
and with the following couplings: nearest-neighbor (NN) 
couplings K x , K y , K z along the cubic symmetry axes, 
next-nearest-neighbor (NNN) couplings Ji,J%, J3 only in 
the ±(1, 1, 0), ±(0, 1, 1), ±(1, 0, 1) directions [but not in 
the ±(-1,1,0), ±(0,-1,1), ±(-1,0,1) directions] , and 
a third-NN coupling K only in the diagonal ±(1,1,1) di- 
rection (Fig. 3). The corresponding anisotropy matrix is 



K, 




FIG. 3: Lattice points Xj of the primitive cell (cube) of the 
anisotropic simple-cubic lattice model (|2.ip and (|2.7p with 
A/l. Solid, dashed, and dotted lines indicate the NN cou- 
plings K a , the NNN couplings Ji, and the third-NN coupling 
K. 

obtained from (|2.10[) as 



A = 2d 



-2 



/ D x _ ,h+K J 3 + K 
J!+K D y _ J 2 + K 
\J 3 + K J 2 +K D z 



(2.14) 



with the diagonal elements D x = K x + J\ + J3 + K , 
D y = Ky + J 2 + Ji+K,D z =K z + J 3 + J 2 +K. For 
quantitative analytical and numerical studies this model 
with seven different couplings would, of course, be much 
too complicated. We shall present explicit quantitative 
results only for two nontrivial cases: 

(i) Model with three-dimensional anisotropy: isotropic 
ferromagnetic NN couplings K x = K y = K z = K > 
and three equal anisotropic NNN couplings J\ = J 2 — 
J3 = J. MC simulations for three-dimensional Ising 
models with this type of anisotropy (with K = 0) 
have been performed by Schulte and Drope [42| and 
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by Sumour et al. [43|. The corresponding reduced 
anisotropy matrix (with K ^ 0) is 



A = (1 — 3w + 2w 6 ) 



3X-1/3 



1 w w 
w 1 w 
w w 1 



which depends only on the single anisotropy parameter 

J + K 



w 



K + 23 



A" 



(2.16) 



The eigenvalues of A and A are Ai = 2a 2 (K + 4J + 3K) , 
A 2 = A 3 - 2_a 2 (K + J) , and Ai = (1 -3w 2 + 2w 3 )~ 1 / 3 (l + 
2w) , A 2 = A 3 = (1 — 3u> 2 + 2u> 3 ) _1 / 3 (l — w) , respectively. 
The eigenvectors 




, e 



(2) 




,0) 



(2.17) 

defining the principal axes are not parallel to the cu- 
bic symmetry axes. The possible range of w consis- 
tent with detA(iu) > is — \ < w < 1. In the limit 
K — > 0, J — > at fixed A ^ corresponding to to — > 1 
the model describes a system of variables ^ on decoupled 
one-dimensional chains with NN interactions K. In the 
previous version [l2j of this model with K = the range 
of w was restricted to — ^ < w < ^. A vanishing of A 2 
and A3 occurs for J — ► —A. At some value w = wlp near 
— 1/2 (corresponding to Ai = 0) our model is predicted 
to have a Lifshitz point with a wave- vector instability in 
the direction of eW, i.e., in the (1, 1, 1) direction (see also 
Sect. VIII. E). 

(ii) Model with two-dimensional anisotropy: An 
anisotropic NNN coupling J\ = J is taken into account 
only in the x — y planes of the three-dimensional sc lat- 
tice whereas all other anisotropic couplings J 2 , J3 and K 
vanish. This model is of interest for comparison with 
the MC data by Selke and Shchur for the two- 
dimensional anisotropic Ising model as will be discussed 
in Sect. VIII. For further recent studies of the anisotropic 
two-dimensional Ising model see also [72| . 

The bulk critical behavior of both models (i) and (ii) 
belongs to the same d = 3 universality class as that of 
the isotropic model with K x = K y = K z = A > and 
J\ = J2 = J3 = K = provided that A a > , a = 1, 2, 3. 



B. Rotation and rescaling: shear transformation 

In order to derive an appropriate representation of 
the anisotropic bulk order-parameter correlation function 
(see Sect. Ill), to develop an appropriate formulation 
of finite-size perturbation theory (see Sect. IV), and to 
treat the anisotropic Hamiltonian H by RG theory (see 
Sect. V) it is necessary to first transform H such that 



the 0(k a kj3) terms of 5K(k) attain an isotropic form. 
This is a shear transformation that consists of a rotation 
and rescaling of lengths in the direction of the principal 
axes [IH . The rotation is provided by the orthogonal ma- 



(2.15) trix U with matrix elements U a £ 




where e^ 1 denote the Cartesian components of the eigen- 
vectors e( Q '. The rescaling is provided by the diagonal 
matrix A = UAU -1 with diagonal elements A Q > 0. In 
k space the transformation is k' = A 1 ^ 2 Uk such that the 
0{k' a k'p) term of 8K is brought into an isotropic form 
with A' = 1, 



SK(k) = SKiV^X- 1 ^') = SK'{k') = J2 k 2+°( kli ) ■ 

a=l 

(2.18) 

In real space the transformed lattice points are = 

1/9 

A ' Uxj. This transformation leaves the scalar prod- 
uct k' • Xj = k • xj invariant. Thereby the volume of the 

primitive cells is changed to v' — (det A) _1 / 2 w. Corre- 
spondingly the total volume of the transformed system is 
V = Nv' = (dct A) _1 / 2 V^ with the characteristic length 
JJ = V a / d . 

Our transformation is defined such that the values of 
the couplings Kij on the transformed lattice as well as 
the temperature variable ro(T) including the values of 
ro c ,ao and t are invariant [see also Eq. (|4.32[) below]. 
This requires us to perform the additional transforma- 
tions ip'j = (det A) 1 / i cpj,u' = (det A) _1 / 2 u and 



ti = (detA) 1 / 4 /*. 



(2.19) 



In terms of the Fourier transform </?'(k') — 



v 'J2'jLi e tk ' Xj (p'j the transformed lattice Hamil- 



tonian reads 

k' 

+ <V'^ £ kVq ^'(kV(pOM(qV(-k' - P' - q')] 



-h'<p'(0) 



(2.20) 



We illustrate this transformation by the example of 
the simple-cubic lattice model shown in Fig. 3. The 
primitive cell with the volume v' = (det A) _1 / 2 a 3 of the 
transformed system is shown in Fig. 4. It has the shape 
of a parallelepiped whose lengths and angles are deter- 
mined such that the transformed second-moment matrix 
A' = 1 is the unity matrix although there are still the 
same NN couplings K a , NNN couplings Ji , and third-NN 
coupling K as in the simple-cubic lattice model of Fig. 
3. 

Working with H' rather than H will be of advantage in 
the context of bulk properties and bulk renormalizations 
in Sect. V. This is not the case for the confined system. 
Although the 0(k' 2 ) term of SK'(k') in (|2~T5j) and (j2~^D|) 
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FIG. 4: Lattice points of the primitive cell (parallelepiped) 
of the transformed lattice model (|2.20[1 and (|2.23|l with the 
volume v' = (AiA2A3) _1 / 2 5 3 . Solid, dashed, and dotted lines 
indicate the NN couplings K a , the NNN couplings Ji, and 
the third-NN coupling K. The couplings are the same as in 
Fig. 3 but A' = 1 (compare Figs. 1 and 2 of [Q). 



looks quite simple, namely k' k' with a trivial anisotropy 
matrix A' = 1, the summations are nontrivial. 

For concreteness consider the simplified example (i) 
with the matrix (|2.15p and the eigenvectors (|2.17|) . While 
the k vectors of the sc lattice have the simple form 



k = 





(2.21) 



with the integer numbers m t = 0,±1,±2, the k' 
vectors are considerably more complicated, 



k' = 



k[ 

k' 2 | = 

^3 



2vr 



^ (mi + m 2 + ms) V2A1 

(m 2 — mi)\/3A2 
y (mi + m 2 - 2m 3 )y/% 

(2.22) 

These k' vectors reflect the shape and lattice structure of 
the transformed system. Thus the price paid for trans- 
forming A ^ 1 to A' = 1 is to work with more com- 
plicated k' vectors. This example demonstrates that the 
effect of anisotropy cannot be eliminated for confined sys- 
tems. In our applications the summations in finite-size 
perturbation theory will be performed in the simpler k 
space whereas bulk integrals (with infinite cutoff) are sim- 
plified in k' space. 

In real space the Hamiltonian H' reads 



H' 



N 
i=l 



N 



(2.23) 



By substituting the transformations defined above one 
easily verifies 

H(r ,h,uo,K it j,v,L) = H'(r ,h', u' ,K it j,v',L'). 

(2.24) 

The measure for the temperature distance from criticality 
^0 — ?*0c = dot is the same for both H and H' . Defining 



the free energy density /' (divided by fc^T) as 

/'(*, ti, L') = -V'~ 1 In Z'(t, ti, L') , (2.25) 



Z'(t,h',L') 



11 f v l\n(2-d)/{2d) 



exp (-H') , (2.26) 



one obtains the exact relations 

Z(t,h,L) = (detA)- nN /W Z'(t,ti,L') , (2.27) 



f(t,h,L) = (detAy 1/2 f'(t,h',L r ) + [n/[2dv)] ln(det A) . 

(2.28) 

The last term is a bulk contribution, i.e., independent 
of L. Furthermore it is independent of t, i.e., a non- 
singular bulk contribution, thus the singular bulk parts of 
f b (t, ti) = /(*, h, 00) and of f(t, ti) = /'(*, ti, 00) as well 
as the total singular parts of f(t, h, L) and of f'(t, ti, L') 
are related by 



f s>b (t, h) = (det A)" 1 / 2 f ( t , h') 



(2.29) 



f s (t, h, L) = (det Ay^fit, ti, L') . (2.30) 

The bulk correlation function of the transformed system 
is 

G' b (x[-^;t,ti) = lim {< >' -(< <p' >') 2 } 

J V'—ioo J 

(2.31) 

where < ... >' denotes the average with the weight ~ exp 
{-H'). It is related to G b by 

G fc (x; t, h) = (det A)- 1 / 2 G^A" 1/2 Ux; t, (det Af^h). 

(2.32) 

The corresponding relation between the Fourier trans- 
forms is 

G 6 (k; t, ti) = G|(A 1/2 Uk; i, (det A) 1/4 h). (2.33) 



(In the arguments of Eqs. (|2 . 2T[) - (|2.33[) we have, 
for simplicity, not indicated explicitly the additional 
transformations of uq = (detA) 1 / 2 ^ and of v' = 
(det A)~ x / 2 v.) In terms of the transformed field <p'(x') = 
(det A) 1 / 4 ^(U _1 A 1/2 x / ) the Hamiltonian (|2~T3|) attains 
the form of the standard isotropic Landau-Ginzburg- 
Wilson Hamiltonian 



H field = H' field = |dV[^^(x') 2 + i(VV) 



+ u (ip ) -hip 1 



(2.34) 



where V'<p' = (dip' /dx[, . . . , dip' /dx' d ) with a trans- 
formed cutoff. 
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III. BULK CRITICAL BEHAVIOR OF 
ANISOTROPIC SYSTEMS 



Before turning to finite-size theory of anisotropic sys- 
tems it is necessary to first discuss the bulk critical be- 
havior of anisotropic systems and its relation to that of 
isotropic systems. We are not aware of such a discus- 
sion in the literature. It is well known that anisotropic 
systems and isotropic systems have the same critical ex- 
ponents (in a limited range of the anisotropy, see, e.g., 
|73j |. and references therein). Within the ip 4 theory this 
is immediately seen from the relation of dimensionally 
regularized bulk integrals (at infinite cutoff) such as 

oo oo 

u J(r Q + k- Ak)- 1 = u' J (r + k' • k')" 1 
k k' 

= - ^- 2)/2 W , (3.1) 



and T >T C [12j . Here we extend this discussion to h 7^ 
and T < T c and provide an appropriate formulation of 
the order-parameter correlation function and of the scat- 
tering intensity in terms of both the eigenvalues X a and 
the reduced eigenvalues A Q of the anisotropic system. We 
also present the derivation of several amplitude combina- 
tions for anisotropic systems in terms of universal scaling 
functions. 

All of the thermodynamic bulk relations given in the 
following subsections A and B remain valid also in the 
presence of subleading long-range interactions of the type 
(|2.9[) . This is not the case, however, for bulk correlation 
functions at T ^ T c , h = and h ^ 0, T = T c in the 
large-distance regime corresponding to the shaded region 
in Fig. 2. 



Two-scale factor universality in isotropic bulk 
systems 



OO OO j oc 

/ = (detA) 1 /^ ee (detA) 1 /^ / *£ , 

k' k Q=1 -oo 

(3.2) 

provided that det A > 0. We see that the A dependence 
is completely absorbed by the coupling u' and that the 
d = 4 pole term ~ e -1 does not depend on the ma- 
trix A. This leads to identical field-theoretic functions 
for anisotropic and isotropic systems (as functions of the 
renormalized couplings v! and u, respectively) and yields 
the same critical exponents and fixed point value vl* = u* 
for anisotropic and isotropic systems (see Sect. V). (The 
d = 2 pole of the integral (|3.1|) which has nothing to do 
with the critical behavior in d > 2 dimensions can be in- 
corporated in the geometric factor [62| which is finite 
in 2 < d < 4 dimensions [sec (|5.2I) below].) 

For this reason not much attention has been paid to 
the role played by anisotropy in bulk critical phenomena. 
This is not justified, however, in the context of the im- 
portant feature of two-scale factor universality 0, Q . Its 
validity has been established by the RG theory only for 
isotropic systems at h = with short-range interactions 
[27l [H, [H, [33| . A brief derivation was also given by 
Privman and Fisher [f| and by Privman et al. ]§[ using 
scaling assumptions at h ^ 0. Their ansatz for the order- 
parameter correlation function, however, is not valid for 
anisotropic systems since they assumed the existence of 
a single bulk correlation length ^ . Recently it has been 
pointed out that two-scale factor universality is absent 
in anisotropic systems jl2j and that anisotropy has an 
effect on several universal bulk amplitude combinations 
[HI but no derivation was given. In particular, two im- 
portant universal amplitude relations derived by Priv- 
man and Fisher @ [Eqs. (f3~10|) and (f3~TTj) below] have 
not been discussed in the context of anisotropic systems. 
Furthermore, the bulk order-parameter correlation func- 
tion of anisotropic systems was discussed only for h = 



First we summarize the bulk critical behavior of sys- 
tems described by the (asymptotically) isotropic lattice 
Hamiltonian H' and the continuum Hamiltonian H'j ield . 
Near T c the bulk free energy density can be decom- 
posed uniquely into singular and non-singular parts as 
m h') = fUt, h') + f^ b (t) where f' nsb {t) has a regu- 
lar t dependence. It is well established that f' s b has the 
asymptotic (small t, small h!) scaling form below d = 4 
dimensions 



f' b (t,h')=A' 1 \t\ dv W ± (A 2 h'\t\ 



-0S\ 



(3.3) 



with the universal scaling function W±(z), 00 < z < 00. 
This function is independent of the cutoff procedure of 
H'fieid- W e use the normalization W+(0) = 1. The two 
amplitudes A[ and A' 2 are nonuniversal. 

Because of isotropy it is justified to define a single 
(second-moment) bulk correlation length 



1 Ex- ^G' b (K';t,h') 
2d Ex' G> b (x<;t,h') 



1/2 



(3.4) 



above and below T c , respectively. In (|3.4[) we have as- 
sumed sufficiently rapidly decaying correlations, i.e., gen- 
eral n for T > T c but n = 1 for T < T c . We assume, in 
the asymptotic region |x'| 3> (v 1 ) 1 ^ ,£' ± 3> (v') 1 ^ and 
for |x'|/£j_ < O(l) and small h', the asymptotic isotropic 
scaling form 

G' b (x';t,h') = Dl|x'|- d+2 -' 7 $±(|x'|/4^2^l*r^) . 

(3.5) 



e ± (t,ti) = e 0+ \trx ± (D> 2 h'\t\ 



-0S\ 



(3.6) 



with universal scaling functions $±(x, y) and X±(y). Wc 
use the normalization A + (0) = 1, thus £' + (t,0) = £,' 0+ t ~ v 
above T c . The length £q + will be needed as a reference 
length in the formulation of renormalized finite-size the- 
ory in Sect. V. The corresponding scaling form of the 
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Fourier transform G' b of (|3.5p is 



where D'^ is the asymptotic (large-x') amplitude of 



G' b (k':t,h') = D[\k'\- 2+ ^ ± (\k'\Z> ± ,D' 2 ti\t\-f 3S ) , (3.7) 



$±(a/, 2/') = 27r (d - 1)/2 r((d - l)/2) -1 ^ ds s 1 "" 



x J d(cos$)(sm$) d - 3 e- lscos *<f>±{s/x',y') . (3. 

-l 



The three amplitudes D[, D 2 , and £ 0+ in Eqs. (|3-5[) - 
(|3.7p are nonuniversal. The basic content of two-scale 
factor universality is that all of these amplitudes are uni- 
versally related to the two thermodynamic amplitudes 
A[ and A' 2 . The relations read 

(£o+) d ^l'i = Qx(d, n) — universal , (3-9) 

A' 2 /D' 2 = P 2 (d,n) = universal , (3.10) 

D l 1 {A' 2 )- 2 {A' 1 )- 1 -~< /{dl ' ) = ^(d,n)=um«eraai. 

(3.11) 

In Ref. the universal constants P 2 and P3 were de- 
noted by Q2 and Q3. In order to conform with Refs. 
@> [74| and to avoid confusion we reserve the notation Q 2 
and Q3 for the different universal constants in Eqs. (|3.13|) 
and (|3.14|) below. For the sake of clarity we present the 
explicit expressions for Qi and Pi in terms of universal 
scaling functions in Appendix A. An equivalent formula- 
tion of Eq. PT^|) is 

lim [f' Sjb {t,0) £' + {t,Q) d ] = Qi(d,n) = universal . 

(3.12) 

The validity of (13.91) a nd (|3.12p has been established by 
the RG theory [iJlHO]- 

Furthermore, the following amplitude ratios 



(r'jr'M'JZU) 



2-n 



Q 2 (d,n) — universal , (3.13) 



D' 00 (Zo+) 2 ~yr'+ = Q 3 (d 7 n)= universal , (3.14) 

have been proposed [74| to be universal. The con- 
stants ,T' C and £' c are defined as follows: 0) = 
T' + t-i for t > 0,x' b (0,h') = T' c \ti\-~*/^ where 
X ' h (t,h>) = -d 2 f b (t,h')/dh' 2 , and 4(0, h') = & = 
g ^-u/m. The length £' h with the amplitude £ is a 
natural reference length of finite-size theory at t = 0, h' ^ 
[see (|6.1ip below]. D'^ is the asymptotic (small k') am- 
plitude of G' b (k';0,0) « i5^|k'|- 2 +''. Alternatively Eq. 

dHUU) can be formulated as D' oo (C 0+ ) 2 ^ /T' + = Q 3 {d, n), 
or, equivalently, 

hm {GUx';0,0)(|x'|/^+)^ 2+,? }(^ + ) d /r' + 

|x I — >oo ^ J 

= Q 3 (d,n) = (D'^/D'J) Q 3 (d,n) = universal (3.15) 



G;(x';0,0) « DJJx 



/i-d+2-77 



The derivation of Eqs. 



(|3~T3|) - (|3~T5l) is sketched in App. A. Again, all of the 
constants on the left-hand sides of (|3.13p - (|3.15p are uni- 
versally related to A[ and A' 2 . 

Below T c at h' = we have, for n = 1, = 
with the universal ratio 

= ^-(0) — universal . (3.16) 

Previously the bulk relations (|3.9p - (|3.16p were 
expected to be universal for all systems within a 
given universality class [8[. Consistency with the 
universality of (|3.13p and (|3.14p was found [74[ for 
two-dimensional (square and triangular) Ising lattices 
and three-dimensional (sc and bcc) Ising lattices with 
isotropic nearest- neighbor interactions (see also [751]). All 
of these systems, however, belong to the subclass of 
asymptotically isotropic systems with an anisotropy ma- 
trix A — cqI or A = 1 and with an isotropic scattering 
intensity [74]]. Also the honeycomb-lattice Ising model 
considered in Q belongs to that subclass, with a con- 
stant Co different from that for the triangular lattice or 
the square lattice. As will be shown in Subsect. B be- 
low, Eqs. ([22]) and (|3~T2j) - (|3~To) must be reformu- 
lated for anisotropic systems with noncubic anisotropy 
at 0(k a kp) whereas p.lOp and (|3.1ip remain valid also 
for anisotropic systems provided that A[ ,A 2 ,D[ and D' 2 
are transformed appropriately. 

It has been shown 

0, Ha EH m 

that the univer- 
sal scaling form (|3.5p . (]3.6|) is not valid in the regime 
r' = |x'| above T c . The same can be shown for 

n = 1 in the regime r' ^S> below T c . Note that this 
regime is part of the asymptotic critical region r' 3> a 
and £± ^> a corresponding to the shaded area in the r'~ l 
- ^^T 1 plane in Fig. 2. In this regime, corrections to scal- 
ing in the sense of Wegner [l9[ are still negligible. One 
must distinguish at least four cases: (i) For the (f 4 lat- 
tice model with short-range interactions, the exponential 
decay above T c depends explicitly on the lattice spac- 
ing a via the exponential correlation length £ e [ill l76j ; 
in Sect. X we shall show that it also depends on the 
bare four-point coupling uq. (ii) For the ip 4 continuum 
theory with a smooth cutoff A, the exponential decay de- 
pends explicitly on A via £ e (A) d, [l(| (see also Sect. 
X). (iii) For the (p continuum theory with a sharp cutoff 
A, a nonuniversal oscillator y p ower-law decay dominates 
the exponential decay [ill. l77l|. (iv) In the presence of 
subleading long-range interactions of the type (|2.9p . the 
power law ~ b\x!\~ d ~ a [3,[IB] dominates the exponential 
short-range behavior. For T > T c this has been shown 
explicitly for the mean spherical model where the asymp- 
totic structure for |x'|/£+ ^> 1 at h = is [l6| 



G' b (x';t,0) 



id— 2 





(£) 







Ma-2 



V 



Ax 



withX>(|x'|/£ + ) - 
ing in the sense of 



(3.17) 

(|x'|/£ + ) -4 . In all cases (i) - (iv), scal- 
, (|3.6p and two-scale factor univer- 



12 



sality are violated in the regime |x'| ^> £' ± (shaded area 
in Fig. 2) because, in addition to the reference length 
£q + , the nonuniversal lengths a, u Q 1 , A -1 , and fe 1 /^- 2 ) 
govern the leading large |x'| behavior. 

B. Absence of two-scale factor universality in 
anisotropic bulk systems 

Now we turn to the anisotropic system. According to 
(|2.19[) , (|2.29p and (|3.3p . the asymptotic scaling form of 
f s ^b is given by (jl.ip with the nonuniversal amplitudes 

Ax = Ai(det A)" 1 / 2 , A 2 = A' 2 (dct A) 1 / 4 . (3.18) 

In order to represent the order-parameter correlation 
function (|2 . 1 1[) in an appropriate asymptotic scaling form 
it is necessary to employ both of the diagonal matrices A 
and A with diagonal elements X a and X a . Using (|2.32|) . 
(f2T33|) . f3~5f and ([3~7f we write G b and G b as 

G b (x;t,h) = Di|A _1/2 Uxr !i+2 -' 7 

x $ ± (|A- 1 /2 Ux |/4,^ 2 / l |r /?5 )(3.19) 

G b {k;t,h) = i?i|A 1/2 Ukr 2+ " 

x 8±(|A 1 / a Uk|&,£> 2 fc|t|-' M ) (3.20) 

with the nonuniversal amplitudes 

L>i = D[(det A) ( - 2+7l)/(2d) , (3.21) 

D 2 = D' 2 (dct A) 1/4 . (3.22) 

Here we identify the spatial variable r' in the scaling 
argument of $±(r'/£±,0) used in Fig. 2 as 

r' = |x'| = |A- 1/2 Ux| , (3.23) 

which, for given x, depends on all of the d(d + l)/2 pa- 
rameters contained in A. 

While the simple transformations (|3 . 18|) and (I3.22[) fol- 
low immediately from the transformations of h,ipi, and 
V, the transformation of D\ is less trivial. Using Eqs. 
pTTg]) . ([3T21]) . and (|3~22)) we find that the universal am- 
plitude relations (|3 . 1 0[) and (|3.11|) of isotropic systems 
remain valid also for anisotropic systems : 

A2/D2 = P%(d, n) = universal , (3.24) 

D 1 A^ 2 A^ 1 ~' r/(d ' ,) = P 3 (d,ri) = universal (3.25) 

with the same universal constants P 2 and P3 as in (|3.10p 
and (|3.1ip . Eq. (|3 . 19[) differs from the representation of 
G b of (12J at h = where, instead of Z?i, an overall am- 
plitude A' G = D[(det A) -1 / 2 was employed. The latter 
representation is inappropriate as A' G is not universally 



related to A\ and A 2 . The relations ([3~T0)) and (|3~24)) 
follow from the sum rule (see App. A) 

x' b (t, ti) = -d 2 ti(t, h')/dh' 2 = „' £ G^'; t, ti) 

= X b(t, h) = -d 2 f b {t, h)/dh 2 = v^ G b (x; t, ft)(3.26) 

X 

Less obvious are the relations (|3.1ip and (|3.25p . Their 
derivation is, in fact, based on an additional assumption 
about the unsubtracted order-parameter correlation func- 
tion (see App. A). The physical significance of (|3.25p is 
that, at criticality, the bulk correlation function and its 
Fourier transform, if expressed in term of A, 

G b (x; 0, 0) = D 1 $±(0, 0) |A" 1/2 Ux|- d + 2 -" , (3.27) 

G 6 (k;0,0) =Di $±(0,0)|A 1/2 Ukr 2+ " , (3.28) 

have an overall amplitude Dx that is universally deter- 
mined by the thermodynamic amplitudes Ax and A 2 of 
the bulk free energy f s b . Unlike in isotropic systems, 
however, the spatial dependence of G b and the k depen- 
dence of G b at criticality are governed by the d reduced 
nonuniversal eigenvalues A Q (with d — 1 independent pa- 
rameters). In addition, the knowledge of d(d — l)/2 
nonuniversal parameter is needed in order to specify the 
orthogonal matrix U, i.e., to specify the directions e^"-* 
of the principal axes relative to the symmetry axes of 
the system. Thus l + {d-l) + d(d- l)/2 = d(d + l)/2 
nonuniversal parameters are needed at T = T c and h = 0, 
and d{d + l)/2 + 1 nonuniversal parameters at finite h. 
These parameters can be measured by elastic-scattering 
experiments at bulk criticality of anisotropic solids. 

Now we discuss the temperature and h dependence of 
G b away from criticality. Along the direction of the 
principal axis a the spatial dependence of (|3.19p is, for 
|£(«)|/4 a ><0(l), 

G b (5cW;t,h) = D 1 (\x^\/X 1 J 2 )- d+2 ~ T ' 

x *± (\Z^\/^\D 2 h\t\-P s ){Z.2$) 

with x( Q ) = x^e^ and, because of (j2~Ti?j) and ([3~B|) . 

£\t,h) = A^&fetdetA) 1 /**) 

= &\t\-"X ± (D 3 h\t\- fiS ). (3.30) 

Thus along the different principal axes (see Fig. 1 (b) 
of Ref. [1J|) there exist d different principal correla- 
tion lengths (t, h) which constitute d different nonuni- 
versal reference lengths with d nonuniversal amplitudes 

£0+ = ^J 2 Co+- Their ratios 

= (Ao/A,) 1 / 2 (3.31) 
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are also nonuniversal. Below T c at h = we have £_ 
£o-'l^l -I/ with the universal ratio for each a 



(a) 



£o-/Co+ = = universal 



(3.32) 



but for a 7^ /3 the ratios f^/Co- = (WA^) 1 / 2 and 
£o-V£o+ = (A Q /A/3) 1/,2 X_(0) are nonuniversal. Because 

of a; = a x nLi A « /2 ^ E q s - (E3 and im p!y 



AxIIC 



(a) _ 



= lim 



a=l 

= Qi(d,n) = universal 



(3.33) 



The susceptibility Xb(t,h) with x&(i, 0) = T+i 7 above 
T c and x&(0,ft) = T c \h\~^^ s have the amplitudes T+ = 
r' + and T c = r^(det A)->/^ s \ Here we have used 

(|2~ig|) and f3"S5]l . From Eq. (|3~3l]|) we have £^(0, fo) = 

with 

= A* /2 (det A)-"l^ s \' c . (3.34) 
Eq. (|3.13[) then implies for each a = 1, d 

(r+/r c ) (^Veff) 2 "" = Qa(d,») = universal ,(3.35) 



and from Eqs. (|3 . 1 5|) and (|3.29p we obtain for each (3 

d 

n 

a=l 



lim | Gfc (x^;0,0)fe^ 



d - 2 +"i n 



Qz{d, n) = universal 



(3.36) 



Qi , Q2 and Q3 = (D00/-D00) Q3 are the same univer- 
sal numbers for both isotropic and anisotropic systems 
within the same (d, n) universality class. 

Similar reformulations of universal amplitude relations 
are necessary for and Rj involving the surface ten- 
sion, Eq. (2.58) of Ref. and the stiffness constant 
(superfiuid density) p s = £%T d , Eqs. (2.17) and (3.54) 
of Ref. [8j, respectively. Corresponding nonuniversal 
anisotropy effects must be taken into account in the for- 
mulation of universal relations involving correction-to- 
scaling amplitudes (Wegner [19j amplitudes) as well as 
of universal dynamic bulk amplitude combinations [78j 
such as R\,R2,R^, and Rr defined in Ref. [§]. 

For completeness, we briefly mention also those uni- 
versal bulk amplitude relations that do not involve the 
correlation length, as listed in Eqs. (2.45) - (2.48), (2.51), 
and (2.52) of Ref. [||. It is straightforward to show 
that, as a consequence of the scaling structure of f' sb , 
Eq. (13. 3p . and of the universality of the scaling func- 
tion W±(z) that these relations remain valid also for 
anisotropic systems, i.e., they are independent of the 
anisotropy parameters A a p. Consider, for example, the 



asymptotic amplitudes A'± and r' ± of the bulk specific 
heat C' b = d 2 f bs /dt 2 = (A' ± /a)\t\- a and of the bulk sus- 
ceptibility x' b = -d 2 f' s ,b/ dh ' 2 = r'il^r 7 01 the isotropic 
system above and below T c at h' = 0, respectively, 
and, correspondingly, C& = d 2 fb, s /dt 2 = (A±/a)\t\~ a 
, Xb — —d 2 fs,b/dh 2 = r±|i|~ 7 of the anisotropic sys- 
tem. (For x'b an d Xb below T c we consider, for simplic- 
ity, only n — 1.) Their amplitude ratios are given by 
A' + /A'_ = A+/A- = W+(0)/W_(0) and by 



E± 



d 2 W + (y)/dy 2 



d 2 W^{y)/dy 2 



(3.37) 



Thus the nonuniversal parameters A a p drop out com- 
pletely. Corresponding statements hold for the amplitude 
combinations denoted by R x , Rg. Ra in Ref. @. 

A Monte Carlo (MC) study [42[ of the anisotropic 
three-dimensional Ising model appeared to be at vari- 
ance with the universality of the bulk susceptibility ra- 
tio (|3.37p . Subsequent MC simulations [H[ of the same 
anisotropic model on larger lattices, however, are consis- 
tent with the universality of (|3.37[) . 

The analysis of anisotropy effects near criticality can 
of course be extended also to the scaling form of other 
bulk correlation functions such as < tp' (x 1 ^ 2 ip' (x'j) 2 >. 
It can also be extended to the case of general n below 
T c where one must distinguish between longitudinal and 
transverse correlations. Furthermore, extensions of this 
analysis should be applied also to critical dynamics |78j | 
and to boundary critical phenomena [j^, H3| • 

In conclusion, all critical exponents and bulk scal- 
ing functions W±(z) ,$±(x',y') with \x'\ < 0(1), and 
X±(y') of anisotropic systems are universal, i.e., they 
are the same as those of isotropic systems in the same 
universality class. However, as far as the bulk corre- 
lation function Gf,(x;i, h) is concerned, the knowledge 
only of the scaling function <&±(x',y') would be empty 
unless one knows how the arguments x',y' of $± are 
related to observable properties. In particular right at 
criticality, the spatial dependence of G&(x;0, 0) is not 
at all contained in the scaling function but only in the 

factor |A~ 1/2 Ux|- d + 2 -'' [see $TT7)]. This requires the 
knowledge of up to d{d + l)/2 + 1 nonuniversal parame- 
ters. As far as the universal bulk amplitude relations are 
concerned, two-scale factor universality (with only two 
independent nonuniversal amplitudes) is valid only for a 
subset of such relations, namely for those that do not in- 
volve the correlation length [such as _|3"TT0|) . (|3"TTT|) . (13~2"4"|) . 
([3T2"5]) . (|337|) and those of Ref. '8} mentioned above]. 
The other relations [such as (|3.33[) . (|3.35|) , (|3.36[) ] pro- 
vide universal relations between quantities depending on 
up to d(d + l)/2 + 1 independent nonuniversal parame- 
ters, thus seven parameters in three dimensions. This is 
the property of multi-parameter universality referred to 
in Table I. 

Furthermore, for anisotropic systems there exist 
nonuniversal anisotropy effects of the large - distance 
regime of G b (x; t, h) (corresponding to the shaded region 
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in Fig. 2) at t ^ 0, ft = and h ^ 0, t = in combi- 
nation with the nonuniversal nonscaling features of the 
isotropic cases (i) -(iv) mentioned at the end of subsec- 
tion A above. 



IV. PERTURBATION APPROACH IN THE 
CENTRAL FINITE-SIZE REGIME 

A. General remarks 

Consider the transformed Hamiltonian H' with a one- 
component order parameter at h! = in a finite geometry 
with a characteristic length L' in the presence of periodic 
boundary conditions. It is expected that, for short-range 
interactions, there exist three different types of finite-size 
critical behavior of f' s (t, L') — f' s-b {t) where f' s b (t) is the 
singular bulk part : (a) the exponential L' dependence 
~ exp(— L'/£' e+ ) f° r l ar g e L'/£' e+ ^> 1 at fixed tempera- 
ture T > T c with £4+ being the exponential bulk correla- 
tion length above T c , (b) the power-law behavior ~ L'~ d 
for large L' at fixed L'/£± ,0 < L'/£± < 0(1), above, at 
and below T c where is the second-moment bulk cor- 
relation length (|3.4|) . (c) the exponential L' dependence 
~ exp(— L'/£' e _) for large L'/£' e _ » 1 at fixed tempera- 
ture T < T c with being the exponential bulk corre- 
lation length below T c . For a description of the cases (a) 
and (c), ordinary perturbation theory with respect to u' 
of the isotropic ip 4 theory is sufficient. For the case (b), a 
separation of the lowest mode and a perturbation treat- 
ment of the higher modes is necessary minmiii. 

For anisotropic systems, the distinction between the 
regimes (a), (b) and (c) remains relevant except that 
there exist no single correlation lengths £ e +,£;+, 
and £ e _. In this and the subsequent sections we treat 
the case (b) on the basis of the lattice Hamiltonian 
(|2.ip for n = 1, ft = and defer the cases (a) and 
(c) to Section X. The case (b) corresponds to the cen- 
tral finite-size region above the dashed lines in Fig. 1. 
For simplicity we assume a cubic shape with volume 
V = L d and a simple-cubic lattice with lattice con- 
stant a. Now the summations ^ k run over N discrete 
vectors k = (k\, &2, • • • , fed) with Cartesian components 
k a = 2-Kma/L, m a — 0, ±1, ±2, • • • , a = 1, 2, • • • , d in 
the range —it /a < k a < it /a. 

The goal is to derive the finite-size scaling form of the 
singular finite-size part f s of the free energy density of 
the anisotropic system (|2.1[) for n = 1 at h = with an 
anisotropy matrix A. We shall show that, for small \t\ 
and large L in the regime (b), f s has the scaling form 



f s (t,L;A) = L- d F(t(L'/£' f/";A) 



(4.1) 



where the scaling argument is expressed in terms of the 
transformed length L' — (det A)~ 1 ^ 2d ' ) L , rather than 
in terms of L, and where £' 0+ is the asymptotic amplitude 
of the bulk correlation length defined in (|3 .4|) on the basis 
of the transformed Hamiltonian H'. Because of (|2.30l) . 



T is identical with the finite-size scaling function of the 
free energy density of the transformed system 



f' s (t,L';A) = L'- d mL'/t' Q+ ) 1/lJ ;A) 



(4.2) 



The advantage of the transformed system is that its bulk 
renormalizations (see Sect. V) are well known from the 
standard isotropic ip 4 field theory. Thus, in order to de- 
rive the scaling function it is most appropriate to de- 
velop perturbation theory first within the transformed 
system with the Hamiltonian H', ([2~2T)]) and (|2~l2"3"|l . with 
v' = (detAJ-Vag* 



B. Perturbation approach 

It is necessary to reformulate in detail the field- 
theoretic perturbation approach of |66| in the context of 
our anisotropic lattice model in order to correctly iden- 
tify the total finite-size part of the free energy density 
f' s including all temperature independent contributions 
oc L'~ d and to identify the new parts of the theory that 
are affected by the anisotropy. The decomposition into 
the lowest mode and higher modes reads (p'j = $' + a'j , 

$' =L'- d <p'(0)=N- 1 J2 ( Pj' ( 4 - 3 ) 

3 

o' 3 =L'- <i ^e lk '- x '^'(k'), (4.4) 

where L' d = (det A)~ 1 / 2 L d . Correspondingly, the lattice 
Hamiltonian H' and the partition function Z' , (|2.26[) . are 
decomposed as H' = H' Q + H'(&, a') , 



H' Q {r ,u' ,L', $' 2 ) = L' d ( -r $' 2 + u' Q $> 4 ) , (4.5) 



JV 



N 



3 = 1 

K 



(4.6) 



i.j=i 



T/d/2 



'exp{-[^ + f'($' 2 )]} , (4.7) 



f'($' 2 ) 



In 



n 



( V ')y d L' d / 2 



do-'(k') 



x exp 



(4.8) 



where 0"'(k') = <^'(k') for k' ^ and H' is expressed in 
terms of o"'(k'). No terms ~ $' a'j and ~ $' 3 er^ appear 
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in (|4.6p because of J^j a 'j = 0- The integration mea- 
sure J d cr'(k') is defined in Eq. (|B.2|) of Appendix B. 
[The corresponding (but different) functional integration 
/ Da of the (cut-off dependent) continuum model was 
not explicitly defined in Eq. (2.11) of 66].] The quantity 
T'($' 2 ) can be interpreted as a constraint free energy, 
with the constraint being that the zero-mode amplitude 
$' is fixed. The quantity exp[— r'(<I>' 2 )] is proportional 
to the order-parameter distribution function of isotropic 
systems [67} which is a physical quantity in its own right. 
Therefore, in contrast to the e expansion approach of 
[32| and [HI, we shall not expand the exponential form 
exp[-f'(<I>' 2 )] but only f'($' 2 ). The advantage of our ap- 
proach has been demonstrated for the specific heat below 
T c in Refs. [HHI. 

Following [65| and [66| we decompose H'(<f>',cr') = 
H[ + H 2 into an unperturbed Gaussian part 



Hi 



N , 

E r 0L 12 



N K 
9 



(4.9) 



and a perturbation part 



H' 2 = v' ifj [6^($' 2 - <)af + Au'^'af + u' Q af] 

(4.10) 

The crucial point is to incorporate the lowest-mode av- 
erage 



M 2 (r ,u' ,L r ) = 



into the parameter 



J d& $' 2 exp(-JJ^) 

3 

OC 

J dV exp(-ff ) 



r OL (r , Uo ,L') = r + 12 Uo M 2 



(4.11) 



(4.12) 



of the unperturbed part H[ and to treat the term 
6mq($' 2 - M^)af of H 2 as a perturbation. The treat- 
ment of the Gaussian fluctuations a' 2 as a perturbation 
is similar in spirit to an earlier perturbation approach 
for Dirichlet boundary conditions [8l| where part of the 
Gaussian fluctuations of the higher modes were included 
in the perturbation part of the Hamiltonian. The pos- 
itivity of r' 0L > for all ro permits us to extend the 
theory to the region below T c . For finite L' , M' 2 and r' 0L 
interpolate smoothly between the mean-field bulk limits 
above and below T r 



lim M' 2 = M' 2 f 





-roHK) 



for r > , 
for r < , 



(4.13) 



lim 



' 0L 



for r > , 



ro ror r ^u, (4W) 
-2ro for ro < . 



The contribution of H[ to L' d r'($' 2 ) is (compare Eq. 
(|B.4|) in Appendix B) 



In 



n 



( v 'y/dL'd/2 



exp(-H[) 



N — 1 
2L ,d 



ln(27r) 



LE ln {[-0L + ^'(k')](«') 2/d } 



2L' d 



(4.15) 



The leading contributions of the perturbation term 
6u' (<f>' 2 - M 2 )af of H' 2 to L'- d f '($' 2 ) read 

6^($' 2 - M"o 2 )5i(ro L ) - 36u 2 (<i>' 2 - M' 2 ) 2 S 2 (r' 0L ) 

+ OK 3 ($' 2 -M^ 2 ) 3 ) (4.16) 



where 

S m {r'^)=L'- d Y J {[r'c 



OL 



SK 



'(k')]} 



(4.17) 



The terms ~ u' Q $>'cr'? and wd "^ 4 of H' 2 yield higher-order 
contributions of O(u 2< £>' 2 ,u' ) which will be neglected 
in the following. We emphasize, however, that leading 
finite-size effects caused by the four-point coupling u' are 
taken into account in Eq. (|4.16|) as it contains the coupling 
between the fluctuations <I>' 2 — Mg 2 of the lowest mode 
and those of the higher modes a'{k.'). For a discussion of 
the order of the neglected terms see also Refs. |66l |81|. 

The starting point for our perturbation expression of 
the bare free energy density (|2.25|) is 



f d/2 °° 



1h(2tt) 



N-l 
2L' d 

2 k^ ln { 



'OL 



5K 



'(k')] {v') 2,d ) 



6^^1(4) - 36< 2 M 4 S 2 (r 0L ) , 



(4.18) 



H 'eff 



-id 



2 



u' eff & 4 



(4.19) 



r' eff =r + Uu'oSiir'^) + U4u' 2 M' 2 S 2 (r' oh ) , (4.20) 



u 'o ff = u o~ 36u 2 S 2 (r' 0L ). 



(4.21) 



Apart from the different form of the lattice interaction 
5K'(k') and the different vectors k', Eq. flTTB"]) differs 
from the previous Eqs. (4.3), (4.11) and (4.12) of the 
isotropic field theory of Ref. [66[ in two respects: (i) In 
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(|4.18p there are additive logarithmic finite-size terms pro- 
portional to L'- d \n[{v') 1 / d ]; in the regime (b) mentioned 
above, they will cancel each other, and a dependence on 
ln^v') 1 /^] will remain only in the bulk part [see l|C.l[) in 
App. C and (|P5)) - P"55l> ]. (ii) In (j4l8]) there are the 
additive logarithmic finite-size terms 



1 

2v 



7 ln(27r) 



1 2tt 
2Ud 17* 



(4.22) 



where v' — L' d /N. These terms are independent of t and 
h, therefore such terms do not affect the physical quan- 
tities considered in Ref. [66| which are derivatives of the 
free energy with respect to t and h. These terms, how- 
ever, must not be omitted in the calculation of the free 
energy itself. While the first term on the r.h.s. of (|4.22|) 
is an unimportant nonsingular bulk part, the second term 
yields a nonnegligible contribution to the universal value 
of the finite-size scaling function T ex at T c which is a 
measurable quantity. (The second term affects the argu- 
ment of the first logarithmic term of the scaling function 
at T c given in (|6.12|) below.) Omission of this term would 
cause a misidentification of the finite-size scaling func- 
tion of the excess free energy density. This would yield 
an incorrect result in a comparison with Monte Carlo 
data [H, HH, H(| that measure the total amplitude of the 
L~ d term of the excess free energy of two- and three- 
dimensional spin models. 

Our approach incorporates, in an approximate form, 
the effect of the finite-size fluctuations $' 2 — Mq 2 of the 
lowest mode amplitude around its average Mq that are 
present in the central finite-size critical region. This is 
not taken into account in the effective Hamiltonian of 
|64| which contains fluctuations of <I>' 2 around zero. Set- 
ting Mq 2 = and r' QL = r in Eqs. P~T5|) - P~2"Tj) would 
yield the bare free energy density corresponding to per- 
turbation theory based on the effective Hamiltonian of 
[64| . This would restrict the theory to the regime tq > 0. 
A foundation of Eqs. (|4.18[) - (14.21)) can also be given 
on the basis of the order-parameter distribution function 



C. Improved perturbation expression 

In its present form the saddle point contribution of the 
lowest-mode integral in (|4. 18)) for large V below T c is 



(«0 

16«i 



'eff 



which, after expansion of (u'q^) 1 with respect to u' , 
would produce arbitrary large powers of u' . On the other 



hand it is clear at the outset that, because of neglecting 
the terms ~ m $' 2 erf and ~ u 'o°~f 01 H' 2l the neglected 
terms in (|4.18p are bulk terms of O(u' ) corresponding 
to two-loop terms. Therefore it is necessary to further 
improve the perturbation expression (|4.18[) . Here our 



reformulation of the In J d& 



, e -H*ft 



term will be guided 



by the requirement that higher-order powers of uq are 
neglected already at the level of H' e ^ , before integrating 
over <&'. For this purpose we rewrite the logarithm of the 
integral over the lowest mode as 

{T Id/2 °? 
J -p [~ H ' eff ] 



In / ds exp 



1 m 

2 



r' eff L' d / 2 
2« // )V2 

L ld/2 



2 4 

s — s 



(w') 2/d (uo e// ) 1/2 



(4.24) 



For the reason given above it is appropriate to expand the 
factors (u'q^)~~ 1 / 2 in both terms of Eq. (|4.24[) in powers 

of Mo and to neglect terms of O(it ) corresponding to 
a truncation of the expansion 

« ff r 1/2 = u - 1/2 + l< /2 S 2 (r 0L ) + 0{uT) ■ 

. (4-25) 

In summary our improved perturbation expression for 
the bare free energy density reads 



i 



2L' d 



^ln{[r OL + ^'(kO] (v') 1/2d } 

k'/O 

oo 

6« 2 5r(r 0L ) 



2L' d 



In 



L fd/2'eff 



v '2/d 



mu' 2 M^S 2 (r' 0L ) 



(4.26) 



with 



'eff 

Vo 



' d ' 2 u'- 1/2 {r [l + 18u' S 2 (r> 0L )] + l2u' S 1 (r' 0L ) 
' 2 M< 2 S 2 (r> 0L )}, (4.27) 



L 

144m, 



'eff 1-1/2 n 



18M 5 2 (r 0L )] 



(4.28) 



( 4 - 23 ) Now, because of u' M(? ~ O(u 1/2 ) at T c , y' eff and 



Jeff 



and the last two terms in (|4.26[) contain terms 
only up to O(u ) at T c . One can verify that in the 
bulk limit below T c the last two terms — 6u A/o 2 S'i and 
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-3Qu' 2 M^S 2 of (|4~26|) which are of 0(1) are exactly can- 
celled by the O(l) terms of the saddle-point contribution 
-(y' eff ) 2 /(16L' d ) of the integral term of g^SJ). Thus 
Eq. (|4.26[) correctly contains the bare bulk free energy 
density = lirriL/^oo /' in one-loop order [i.e., up to 
0(1)] 

C = - ^ + \ J m{[r +^'(k')]K) 2/d }+OK), 

k' 

(4.29) 

fiT = \^M'l s + u' M£ f - ^ 
+ ln{[-2r + 5K'(k')}(v' f /d } + O(u' ) (4.30) 



above and below T c , respectively, where 



/ = (detA)W ^ (detA)^f[ 7 *2 



-7r/a 



(4.31) 



We shall rewrite (|4.26[) in terms of ro — r$ c where 
1 f 1 



roc 



12Un 



iSX'fk' 



12u 



5iT(k) 



(4.32) 



is the critical value of ro up to 0(uq). On the level of 
bare perturbation theory, the application of Eq. (|4.26p 
will be in the central finite-size regime \ro — rod ^ 
0(u'^ 2 L'~ d l 2 ). On the level of the asymptotic renor- 
malized theory this will correspond to the finite-size 
regime < \t{L' /£ ) 1 / v \ < 0(1) above, at, and below 
T c , i.e., the regime (b) mentioned above. If applied to 
the regime L'/£'_ ^ 1 below T c , /' also contains bulk 
and finite-size terms of O(u' ) which would need to be 
complemented by two-loop calculations. 

The right-hand side of Eq. (|4.26[) can be decomposed 

as 

f'(r - r 0c , u' , L', K id ,v') = f^ b (ro - r 0c , u' , K hj ,v') 
+ Sf'(r -r 0c ,u' ,L',K itj ,v'), (4.33) 

where f„g \ is a non-singular bulk part up to linear order 
in r - r 0c , 



f'nsl( r o -r c,u' ,K id ,v') 



ln(27r) 
2v' 



+ \ ln{[5K\V)]{v'Y /d }+ r -^P± 



[SK'(k> 



(4.34) 



As expected from bulk theory |82| , the remaining finite- 
size part Sf has a finite limit for v' — > at fixed tq — r$ c 



in 2 < d < 4 dimensions. It turns out that the resulting 
function depends only on A rather than on A, 

lim 8f'{ro-ro c ,n' ,L',Kij,v') = Sf'(r Q -r 0c , u' , A) . 

(4.35) 

The r.h.s. of (|4.35p can be further decomposed as 

Sf'(r - r 0c , u' , L', A) = /„ s 2 J,(r - r 0c , u' ) 

+ f s (r -r 0c ,u' ,L',A) , (4.36) 

where is a non-singular bulk part proportional to 

{tq — ?*o c ) 2 [82]. We are interested in the asymptotic sin- 
gular finite-size part f s . In the limit v' — > our result for 
Sf does not contain an L dependent non-singular part. 
The limit v' — ► is justified in the power-law regime (b) 
mentioned above where the v' dependent terms of our 
perturbation expression (14.26|) give rise only to correc- 
tions to scaling. However, although the limit (|4.35p does 
exist in the exponential regimes (a) and (c), it is not jus- 
tified to neglect the v' dependencies in the exponential 
arguments, as will be discussed in Section X. 



D. Bare perturbation result 

The calculation of Sf is outlined in App. B and C 
for the power-law regime |ro — rod < O(u' L'~ d / 2 ) , L' ^> 
(v') l ' d ,|r - r 0c | 1/2 (u') 1/d < 1. The result reads for 
2 < d < 4 



Sf(r - r 0c , u' , L', A) = 



-A d (r' 0L )~ e/2 



W L ) 2 , (r -r 0c ) 2 



U 



Ae 



- 18u' 2 M^ 



±1 - \n [ ds exp(-iy^(A) S 2 - S 4 ) 



In 



'eft i 
(i') £/2 



2irw' eff (A) 



+ 2<W L L' 2 ,A) 



3u' M?L> 2 

^2 ll{ r 0L L i A ) 



9u ,2 MfL' 4 
4^4 



l2(r' 0L L' 2 ,A 



(4.37) 



with 



y' eff (A) 



V 



1/2 



{ro - r 0c ) 



-e/a 



l( A d {d-2) r , 
18u °l^^ (r ° L) 1 16^ 

+ 12<(- ^(r 0L )(^ 2 )/ 2 + ^g^/, / \ 



r / 2 (r 0L L' 2 ,A 



14«(^^(^)-/ 2 



16tt 4 



7 2 (r 0L L' 2 ,A) 



(4 
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(A) = 



1-1/2 



1 + 18it 



A d {d-2) 
2e 



16tt 



■/2(4i' 2 ,a; 



r OL( r o - r 0c , u' , L') =r - r 0c + 12u' q Mq, 



m? = (Vu' )-^ Mv'o), 

2/ = (r -r 0c )(y'K) 1/2 , 



J (r 0L L' 2 ,A) 



J ds s r 
o 



exp(- 



2V0 S 



J ds exp(-±y s 2 
o 



dy 
y 



exp 



'0L J 



4tt 2 



-e/2 

(4.39) 

(4.40) 
(4.41) 
(4.42) 

(4.43) 



{(n/y)^ 2 - K d (y, A) + l} - exp(-y) 



V. MINIMAL RENORMALIZATION AT FIXED 
DIMENSION 

The bare perturbation form of 5f requires additive 
and multiplicative renormalizations, followed by a map- 
ping of the renormalized free energy Sf' R from the criti- 
cal to the noncritical region where perturbation theory is 
applicable. It is well known that, for the multiplicative 
renormalizations, the usual bulk Z factors are sufficient 
[8~4| . For both multiplicative and additive renormaliza- 
tions, the absence of L-dependent pole terms has been 
checked explicitly up to Oju'r, 2 ) for the case of periodic 
boundary conditions [HI, [8l| . In particular, there is no 
need for an L-dependent shift of the temperature vari- 
able. We employ the minimal subtraction scheme at fixed 
dimension 2 < d < 4 without using the e expansion [6^ |. 
This approach has already been successfully employed in 
previous finite-size studies [65|, [66|, H3, [H, [6^, and is 
applicable above, at, and below T c with the same renor- 
malization constants. Thus it permits us to derive a sin- 
gle finite-size scaling function of the free energy in the 
central finite-size critical region above, at, and below T c . 

The multiplicatively renormalized quantities are 



(4.44) 



A,Z 



(5.1) 



I m (r' QL L ,2 ,A) = J dyy m -^M-r^L l2 y/{^)] 
o 

x {K d (y,A) — {■n/y) d ^ 2 — 1} , 

(4.45) 

K d (y, A)=J2 exp(-yn • An) (4.46) 

n 

with n = (m, ri2, TI4) ,n a — 0, ±1, ±oo. The be- 
havior of the functions Jo , Ii , and I2 for small and large 
arguments r 0L L' 2 is given in App. C. 

The crucial information on the anisotropy is contained 
in the sum (|4.46[) . By means of the Poisson identity [83[ 
[see also (|B.9[) ] one can show that this function satisfies 

K d (y;A) - (det A)" 1 / 2 (£j ' K d (^-,A- X 

(4.47) 

The sum (|4.46p could formally be rewritten in k' space as 
exp(— y'k' -k') with y' = yL' 2 /(4TT 2 ) but in practice 
there is no advantage of using the more complicated k' 
vectors (see the example (|2.22p in Sect. II). For this 
reason the sums in the three-dimensional calculations in 
Sect. VIII will be performed in k space. 

As expected, the bare perturbation result (|4.37|) for 5f 
does not yet correctly describe the critical behavior: (i) 
In the bulk limit at t ^ the small-i behavior is 8fl ~ 
\t\ d / 2 rather than ~ \t\ du . (ii) At t = the leading large- 
V behavior is Sf ~ L'- d2 / 4 rather than - L'- d . These 
defects will be removed by turning to the renormalized 
theory. 



= Z , 2 ip' with an 



and r = Z r (ro — r$ c ) = at 

arbitrary inverse reference length //. L' is not renormal- 
ized. Furthermore, the reduced anisotropy matrix A is 
not renormalized either as it does not change the ultra- 
violet behavior at d = 4. If our calculation is extended 
to a finite external field h', (|2.19[) . the additional renor- 



malization h' T 



is necessary [6£ 



The geometric factor of bulk theory [see (|3.ip ] 



A d = 



r(3-d/2) 

2 d-2 nd /2( d _ 2 ) 



S d T(l + -)T(l 



(5.2) 



appears naturally in Eqs. (|4.37|) - (|4.39|) rather than the 
more commonly used factor S d = 2 1 ~ d n~ d / 2 {T(d/2)]~ 1 . 
The perturbation results of amplitudes and scaling func- 
tions dep end on the choice of the geometric factor in 
(|5.ip [62[ HI, [66|, [H, [85[ (see, e.g., the universal ratio Qi 
in (|6.19p below, see also the comment after f|5 . 1 6[) below). 
The advantage of the factor (|5.2p is that it describes the 
full d dependence of single-loop integrals in 2 < d < 4 
dimensions such as (|3.ip . in contrast to the factor S d . 
For this reason we have incorporated A d in the defini- 
tion of the renormalized coupling u 1 , (|5.ip . Any other 
choice, such as S d instead of A dl would introduce arti- 
ficial d dependencies into the perturbation results. For 
the same reason we employ A d in the definition of the 
multiplicatively and additively renormalized free energy 
density 

f R (r,u',L',»,A) = 5f'{Z r r,if Z^Z-fA^u' ,L' ,A) 



yT e r 2 A d A(v! ,e) 



(5.3) 
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Because of relations such as (I3.1|) the Z factors Z r (u', e), 
Z u i (u', e), and Z v r(vl ', e) depend on u' in the same way as 
the usual Z factors depend on u in the standard isotropic 
ip 4 theory. The same statement holds for the additive 
renormalization constant A(u' \e) because of 

oo oo 

J ln(r + k • Ak) = (det A)" 1 / 2 J ln(r + k' • k') 

k k' 

= -(detA)- 1 /^^. (5.4) 

Thus the renormalization constants read up to one-loop 
order Z r (u',e) = 1 + 12u'/e, Z u ,(u',e) = 1 + 36u'/e, 
Z v i(u\e) = 1, A(u',e) = -2/s. The Z factors Z u > and 
Z r are sufficient to renormalize y^* , (|4.27p , and w'q^, 
(j4.28|) , whereas the additive renormalization constant 
A(u',e) is needed to absorb the pole term ~ — Ad(ro — 
r 0c) 2 /(4e) in the square brackets of (|4.37p . After substi- 
tuting these renormalization constants one verifies that 
the resulting renormalized free energy density f' R has a 
finite limit for e — > at fixed u' > 0. 

We define the dimensionless amplitude function 



F R (r/v 2 ,u\ L'fi, A) = n~ d Affair, u', L', A) 



(5.5) 



From the \x independence of 8f'(ro — ro c , u' , L', A) one 
can derive the renormalization-group equation (RGE) for 
the amplitude function 



(p d„ + r( r d r + f3 u> d u . + d)F' R {r/ii 2 ,u', L'p, A) 



[r 2 /(2p i )]B{u') 



(5.6) 



Eqs. (|4.37p - (|4.39|) the parameter r' 0L will appear in the 
form of the effective renormalized counterpart 

r' h (l) = r'M^^A^u'^L') 

= r(0 + 12(^A- 1/2 u'(l)^(Lr d/2 My'(l)) , 

(5.8) 

y'(l) = r{l)p-H-\L'M d /*A)!\'(l)-W. (5.9) 

Correspondingly, the effective renormalized counterparts 
of y l Q and of w'q appearing in the renormalized form 
of the logarithmic part of Sf are given by 

y>eff(l,A) = (l»Ly/*A 1 J 2 u>(l)-^ 




l + 18u'(l)R 2 ( 
r' h {l) 



1(0 
p 2 l 2 

l(iL',A 



lfiL',A 



xR 



p 2 l 2 

iu(^LT d/2 A- 1/2 u'(i) 3 ^My'(i)) 



n 2 i 2 



,lfj,L',A) 



(5.10) 



■j' eff (l,A) = u'(l)-^ 2 \l + 18u'(l)R. 



p 2 l 2 



, IpL', A 



(5.1 



where the field-theoretic functions f3 u i(u\ e), Cr(V), and 
B(v!) are defined as usual [H, Eq. (|5.6[) . however, 
differs from the corresponding bulk RGE (119) of [HJ 
since here we are using r rather than the bulk correlation 
lengths £,'± as the appropriate measure of the temperature 
variable. Using r rather than £± is of advantage in finite- 
size theories where a single finite-size scaling function 
is derived for both r > and r < 0. The functions 
( r (u') , (3 u i{u',l) and B(u') as well as the fixed point 
value u'* — u* are accurately known (63l . |86| from Borel 
resummations. Integration of the RGE yields 

n(^.»',iv,A) = c |n(^,« , (i),M«,A) 

+ ^/ B ("'<M^/[ 2 «"'<n)-#}y}. 

1 I ' 

(5.7) 

with an as yet arbitrary flow parameter I and u'(l) = 
u' . The effective parameters r{l) and u'(l) are defined as 
usual |62|. 

Eqs. (|4.40p and (|5.7|) show that in the arguments of 
the functions Jo,Ii,Is and of the pole terms ~ e~ x of 



Ri(q,p,A) = e~ 1 q[l-q~ e ^ 2 ]+p e ~ 2 (4:Tr 2 A d )~ 1 Ii(q p 2 , A) , 

(5.12) 



R 2 (q,p,A) = - e -i[i- ? -«/2]_i r «/a 

+ p^IGtt 4 ^)- 1 ^^^^) . (5.13) 

This suggests that the most natural choice of the flow 
parameter / is made by 



(5.14) 



It ensures the standard choice in the bulk limit both 
above and below T c [63| 



lim n 2 l 2 = 



L— *oo 



p 2 l 2 + =r{l+) for T>T C 
p 2 l 2 _ = -2r(Z_) for T <T C 



(5.15) 



and appropriately implies pi oc L'^ 1 for large finite L' at 
T = T c . As a natural choice for the reference length pT l 
we take pT 1 = £' 0+ where 



e 0+ = 



Z r (u', £ )a - 1 Q*exp( /" Cr( "'*l ■■ Cr i"" ) du" 



Pu'(u",e) 



(5.16) 
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is the asymptotic amplitude of the second-moment bulk 
correlation length of the isotropic system above T c , as 
defined in (|3.6p . The dimensionless amplitude Q* — 
1 + 0{u* 2 ) = 1 + 0(u'* 2 ) is the fixed point value of the 
amplitude function Q(l, u', d) of the second- moment bulk 
correlation length above T c [63(. Owing to the choice of 
the factor Ad, (|5.2|) , the O(u') term of Q(l,u',d) and 
the 0(u*) term of Q* vanish [63j, [82j, l87|, similar to the 
vanishing of the order-parameter amplitude function at 
0(u') [63]. The same observation was recently made for 
the correlation-length amplitude within the e expansion 
[5l| where the same geometric factor Ad was employed 
(apart from a harmless factor of 2). In three dimensions 
the amplitude Q* it is accurately known from Borel re- 
summations [13] ■ 

Eqs. (|5.8[) and (|5.14[) determine I = l(t,L') as a func- 
tion of the reduced temperature t and the size L' . With 
this choice of Z, Eqs. (|5.7[) and (|5.5[) provide a mapping of 
the functions F' R and from the critical to the noncrit- 
ical region. 

In summary, the singular part of the contribution Sf 
to the free energy density of the isotropic system is con- 
tained in 

f R (r,u',L',fi,A) = 6f(Z r r,fi e Z u ,Z-^A^u',L',A) 
-^r 2 A d A(u',s) = f R (r(l), u'(l), l»,L', A) 



above and below T c , respectively, where Z+ and Z_ are de- 
termined by Eq. (|5.15p . The last integral term in (|5.17p 
contains both a contribution oc t 2 l~ a / v to the singular 
finite-size part f s and a contribution oc t 2 to the nonsin- 

gular bulk part f'J£\ of Sf [see Eq. (|4.36p . compare also 
Eqs. iJBISJl and l|OU|) ]. 



VI. FINITE-SIZE SCALING FUNCTION OF 
THE FREE ENERGY DENSITY 

A. Result in 2 < d < 4 dimensions 

In order to derive the finite-size scaling function T we 
consider (f5T7|) - ([5T20]) the limit of small I < 1 or I — > 
0. In this limit we have u'(Z) — > it'(0) = u'* = u* , 



(Z) -> y = xQ* (nlL')- a /^A d /2 u*- 1/2 , (6.1) 



5 = t(^') V,y = KL'/^,) 1 /" ■ (6-2) 
In (|6.ip we have used the hyperscaling relation 

2 - a = dv . (6.3) 



A d r(l) 2 
2(l/x)= 



5 



(u'(i')){ 



exp 



2Cr(«'(i")) 



dZ" -i dz' 
(5.17) 



Because of the choice (|5.14p . Eq. (|5.8p implies ^ZL' 
Z = l(x) where l(x) is determined implicitly by 



f R (r{l),u'{l),ln,L',A) = 
- A d (lfJ,) d /(4d) + l8u'{l)L'- d 

s OO 

+ J7d\ - ln / dsexp[-- 



y 



l 



in 



— In / ds exp [ - 

— OO 

2irA 1 d /2 w /e ^(l,A) 



(i^uy/ 2 



iMy'm 2 

y' eff (l,A)(l)s 2 - S 
+ l -Ml 2 H 2 L' 2 ,A) 



= x Q* r«/P")Aj /2 u*- 1/a 



(6.4) 



(6.5) 



Simultaneously, these two equations determine y = y(x) 
Furthermore we have 



■/ eff (l,A)^ W{x,A) = 



u 



-1/2 



l + 18u*i? 2 (l,/,A) 
(6.6) 



3(«^L') £/2 «'(0 1/2 



2n 2 A 



1/2 



tf 2 (2/(0)/iGV£' 2 ,A) 



y' e f f (l, A) -> F(x,A) = f d / 2 4 /2 M *- 1/2 

xJ g'JT 1 ^ [l + 18ti*R 2 (lJ,A)l + 12u*Ri(l,Z~ A) 



9(lnL') e u'(l) 

47T 4 A d 



[M y' m 2 hH 2 , 2 L'\A) . (5.18) + 144u *3 /2[ _ d/2j4 - 1/2i?2( - )ii2(1)U) 



(6.7) 



In the bulk limit the function f Rb (l±/j,,u' (Z±)) = 
liniL^oo f' R (r(l), u'(l), Ifx, L\ A) becomes independent of 

A 

- A d (Z +M )<7(4d) , (5.19) 



The asymptotic (Z — > 0) behavior of the integral in (|5.17p 



f+ b (l +l i,u\l + )) 
f'- b (l^,u'(l^) 
-Ad(l-ti) d 



B(«'(0){exp / [2C r (u'(Z"))- £ ]|^}f 



64u'(Z. 



1 81 //, \ 
Td + 64 u (L) 



- 73(u*) + 0{l a/v ) 
a 



(6.8) 



(5.20) 



is known from bulk theory [8^]. In (|6.8p the subleading 
term 0(l a ^), together with the prefactor r(l) 2 /(lfi) E in 
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(|5.17[) , contributes to the regular bulk term f'^ 2 \ propor- 
tional to t 2 of Eq. (036]). 

In summary, the asymptotic form of the singular part 
f' s l|4.2p of the reduced free energy density of the isotropic 
system at h! = is obtained from f' R , (|5.17p , in the limit 
of small I as 

= J/- d .F(z, A) (6.9) 
where the finite-size scaling function is given by 
¥ vQ* 2 x 2 t a ^ 



F{x,A) = -A d 



18m* [0 2 (y)} 2 



-id 



2a 

'27rAyV(^,A) K 



oo 

In J ds exp [ — ^Y{x, A)s 2 — s 4 ] 



1 op/2 * 1/2 

2 27r 2 Ay 



9 l e 



47T 4 A d 



[tf 2 (y)] 2 / 2 (P,A) 



(6.10) 



This result is valid for 2 < d < 4 in the range L' ^S> 
a and < \x\ < 0(1) above, at, and below T c (but 
not for the exponential regime |x| 3> 1, see Sect. X). It 
incorporates the correct bulk critical exponents a and v 
and the complete bulk function B(u*) (not only in one- 
loop order). There is only one adjustable parameter that 
is contained in the nonuniversal bulk amplitude £q + of 
the scaling variable x, (|6.2|) . For finite L',f' s (t,L') is 
an analytic function of t near t = 0, in agreement with 
general analyticity requirements. From previous studies 
at finite external field [69l . [88j we infer that the extension 
of (|6.9|) to hi 7^ has the structure 

f s (t, ti, L>) = L'- d T(x, ti(L'/0^;A) (6.11) 



where £' c is defined after ()3.14j) . Thus the constants C[ 
and C' 2 in (|1.3p and (|1.4p can be chosen most naturally 

Of particular interest is the finite-size amplitude 
.F(0,A) = T C {A) at T c , 



•^c(A) 



18 -§)u* [^ 2 (0)] 2 -iln 



'27rAy 2 W c (A) 



f c 



oo 

In /" ds exp [ - ir c (A)s 2 - s 4 ] 



1 / 2 



2 ^(^,A)-^/ 1 (/ 2 ,A)-^ I / 2 (f 2 ,A) 



y c (A) = 144u*l? a (0)|Jii(l,Zc,A)+ R 2 (l,l c ,A)\ 

(6.14) 

withtf 2 (0) = r(3/4)/r(l/4) and 

Ri(l,! c ,A) = l 2 ~ d (4ir 2 Ad)^ 1 Ii(l 2 , A) , (6.15) 



i? 2 (l,f c ,A) = -- + ^(167r 4 A d )- 1 / 2 (f 2 ,A) . (6.16) 



In the bulk (large limit Eqs. ([p3)l and |Of])) yield 
lim L ^ 00 f' s (t,L') = f s ± b (t) where 



/r 6 w = 



i 

id 



64u* 



1 81 , 

Ad + 64 U 



2^ *<«*> 



£f—d>dv 



(6.17) 



(6.18) 



above and below T c , respectively, with the universal bulk 
amplitude ratios 



f&mt =Qi 



i 

4d 



(6.19) 



.4- 

3h 



dv l/(64u*) +l/(4d) +81w*/64 + i/B(«*)/(8a) 



l/(4d) + z/S(u*)/(2a) 



(6.20) 



(For Qi compare (|3.12p .) Here we have used the bulk 
identifications 

, f ^+ = 0*^0+^ for T>T C , 
1 ^- = Q*" €o+ a (2|*ir for T<T C , 



(6.12) 



where ~fj 2 = 12u* 1/ ^A^ i/ ^ 2 (0) and 



V2 4-1/2, 



as implied by the choice 15.1 5|) . As noted in Sect. Ill, 
a complete two-loop calculation would yield further bulk 
contributions of 0(u*) in /'^. Owing to the truncation 

(|4.25p . no terms of 0(u* 2 ) and higher order appear in 
(|6~T8j) and (jOO)) . 

In order to present the scaling function 

.F^A) =F(x;A) -Jf(x). (6.22) 

of the excess free energy density f' s ex (t, L'; A) = 
f' s (t,L';A) — f' sb {t) we shall also need the A indepen- 
dent bulk part 



W C (A) = u*~ 1/2 1 + 18 u*R 2 (l, l c , A) , (6.13) 



£ /d /# = Qi^" forT>T c , 
L ld f'- b = Qi\x\ d » for T<T C , 



(6.23) 
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with Qi = (A~ /A + )Qi, representing the large |x| be- 
havior of !F{x,A.). It should be noted that it is not ob- 
vious how to interpret the d x d matrix A for the case of 
non-integer dimensions d. 

In the spirit of the fixed - d minimal subtraction ap- 
proach [62| we shall evaluate Fix, A) and ^-" ea; (i,A) in 
d = 3 dimensions without any further expansion with re- 
spect to u* . This is in contrast to the e expansion which 
is a double expansion with respect to u* and e = 4 — d. 



B. Epsilon expansion 

Considering u* as a smallness parameter and using the 
results of App. C we obtain from (|6.12[) at fixed 2 < d < 4 

2 \ 24^/^(1/4)]^ j 



hi 
1 

"8^ 



rl 



r(i/4) 



12 r(i/4) 
4 /2 r(3/4) 



l 



MO, A) 



4/d 



h{0,A)u* 2/d + 0(u* ,u* i/d ) . 

(6.24) 



Substituting u* — e/36 + 0(e 2 ) and expanding all d 
dependent quantities with respect to e = 4 — d yields the 
e-expansion result at T c 



T C {A) = X - lne + / (A) + A (A) e 1 ' 2 + 0(e) , 



(6.25) 



/o(A) = - - In 18 - In 



2 r (V4) 



1 fdy 



2 J y 

o 



, (6.26) 



A nontrivial question arises if the e-expansion result 
is applied to three-dimensional anisotropic systems with 
a matrix A / 1. It appears that, to some extent, it 
is ambiguous how the physical 3x3 matrix A (which, 
in general, has 5 independent nonuniversal matrix ele- 
ments) should be continued to d = 4 in order to evalu- 
ate the coefficients /o(A) and /i(A). This matrix A in 
(I6.25P is necessarily a 4 x 4 matrix which, in general, has 9 
independent nonuniversal matrix elements, i.e., four ad- 
ditional nonuniversal parameters. It is not unique how 
to choose the magnitude of these four additional matrix 
elements. The results for /o(A) and /i(A) in four di- 
mensions will significantly depend on this choice. 

As a possible choice we propose the following. In order 
to describe the physical system with the symmetric three- 
dimensional matrix 



(6.28) 



with dct A = 1 it seems reasonable to extend this matrix 
to the four-dimensional counterpart 

/ a b c \ 

b d e 

c e f 

\ 1 J 




A = 



(6.29) 



This choice guarantees that no arbitrary anisotropy is 
introduced in the fourth dimension and that det A = 1. 

A corresponding problem would arise in an e expansion 
in d = 2 + e dimensions. In Sect. VIII. D below we shall 
present an example where we compare the anisotropy 
effects of two- and three-dimensional models. The result 
of this comparison supports the suggestion given above 
for the dimensional extension of the matrix A in an e 
expansion. 



C. Large - n limit 



A (A) = 



r(i/4) 

7rr(3/4)V2 



dy 



Ki(y,A)- 



(6.27) 



where now A denotes a 4 x 4 matrix. The e expansion 
result (|6.25|) - (|6.27|) is independent of which renormaliza- 
tion scheme and which kind of perturbation approach is 
used. The same result is obtained if one starts with the 
effective Hamiltonian of Brezin and Zinn-Justin [64l| or 
with the cumulant expansion of Rudnick et al. [321 ] . Be- 
cause of the strict expansion with respect to u* and e, the 
exponential structure of the distribution ~ exp[— H' e ^\ 
is destroyed. As expected, the e - expansion term ~ lne 
is not well behaved for e — > since at d = 4 the finite 
lattice constant a must not be neglected. 



For comparison with the case n — 1 we also present 
the exact result for the finite-size scaling function of 
the free energy density per component /oo(^, L) of the if 4 
lattice model (|2.1I) with v = a d and V — L d in the limit 
n — > oo at fixed uqu. From Eqs. (45) and (46) of Rcf. 
18911 we have 



= /o 



foo(t,L) 



lim {-(nV)~ 1 lnZ(t,0,L)} 

n — >oo 

2^E ln {[x^+^(k)]a 2 X6.30) 



where Z(t, 0, L) is defined by (|2.3|) and \oo i s determined 
implicitly by X ^t = r + 4u ^^" 1 E k [x^ 1 + ^(k)]" 1 . 
The additive constant in (l6~30]) is f = -[ln(27r)]/(2a d ). 
Using the results of App. C leads to the singular part of 
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foe in the regime L' ^> a and < \x\ < O(l) at h 
for 2 < d < 4 



/«,,,(«,£; A) = L- d Too{x;A 



(6.31) 



^o(P 2 ;A) 



-4, 



2(4 - d) 



£P 2 



2 pd 
d 



P 



d-2 



4-d 

~a7 



(6.32) 
(6.33) 




exp 



P^y 
4tt 2 



^(y,A) .(6.34) 



Here 5 = t(L'/^ + )^ v with i/ = (d - 2)" 1 , L' = 
(detA)- 1 /^)^, #, + = (4^7iA d a Q 1 /e) ,y , and u' = 
(det A) _1 / 2 uq. We note that the geometric factor Ad, 



(|5.2j) . appears in (|6.31| - (|6.34j) in a natural way. The 
reason is that only diagrammatic contributions of single- 
loop structures contribute to the large - n limit. The 
function (l/2)Qo(x; A) with x — r^L' 2 is the scaling func- 
tion of the excess free energy density of the Gaussian 
model (see (|B~13| - flEE} of App. B). For T > T c the 
function P(i; A) determines the finite-size scaling form 
of the susceptibility per component in the limit n — > oo 

m 



xto(t,L;A) = L'^g(x;A) }1 /v 



(6.35) 



where g(x; A) = [P(x; A)] 2 . Below we shall present the 
relative anisotropy effect 



Ax+ (A) 



5 (0; A) - ff (0;l) 



(6.36) 



5(0; 1) 

on xto a ^ m three dimensions. 

The result (|6.3ip - (|6.35|) is the extension of the result 
for the isotropic case (see Eqs. (17)-(19) of [ll|) and cor- 
rects Eq. (44) of [HJ where the term — (ln2)/2 should be 
dropped. The scaling function of the excess free energy 
density above, at, and below T c is given by 

T^(x-A) = T x (x\A) - T^ b (x) (6.37) 

with the bulk part 




where Y = (d - 2)A d /[2d(A - 
amplitude is given by 

^• oo (0;A) = ia (P c 2 ;A 



for T > T c , 
for T<T C , 



(6.38) 



d)]. At T c the finite-size 



A, 



where P C (A) 



P 



d-2 



d{4 - d) 
P(0; A) is determined by 
A-d 



P" 



A,, 



(6.39) 



(6.40) 



VII. OTHER FINITE-SIZE SCALING 
FUNCTIONS 

The calculations of the preceding sections can be ex- 
tended to other finite-size quantities. Here we consider 
only those quantities that have been studied in MC sim- 
ulations of anisotropic Ising models [H . l43l . [lij ] . Within 
our </3 4 lattice model (|2.1|) for n = lat/i = 0ona simple- 
cubic lattice with volume V — L d we shall consider the 
susceptibilities x + = V < * 2 >, X - = V(< <I> 2 > - < 
|$| > 2 ), and the Binder cumulant 



i7=l-i<$ 4 >/<$ 2 > 



(7.1) 



where $ = N^Y^jPj ( see j e -g-> @)- These quanti- 
ties remain invariant under the transformation defined 
in Sect. II jl3j |. x ± = (x^Y ,U = V. As a consequence 
we find that, in the regime (b) defined in Section IV. A, 
the finite-size scaling forms of these quantities are 



X 



± (t,L;A) = (**)'(*, L'; A) = (L'/£' )^P±(:r; A) 



(7.2) 



U{t, L; A) = U'{t, L': A) = U(x; A) , (7.3) 

where the scaling functions P5 and U are obtained from 
those of [66| by the replacements Y — > F(ic;A) and 
i?2 — > J?2(lj^A). Note that the functions P^T are 
nonuniversal even for A = 1 since they still contain 
nonuniversal overall amplitudes c proportional to the 
bulk amplitudes of x ± (see also 68]). Here we consider 
only the relative anisotropy effect 



Ax^A) 



P±(0;A) 



P±(0;1) 



^(0;i) 



on the susceptibilities x ± (0, ^;A) at T = T c , 
lytic expressions are in 2 < d < 4 dimensions 



(7.4) 
The ana- 



P+(0;A) = c+ 1 - 18u*P 2 (l, lc, A) 



P"(0;A) = c 



MYc(A)) , 
(7.5) 

-l 



1 - 18u*P 2 (Mc,A) 

{MYc(A)) [MYc(A))} 2 } 



(7.6) 

where the constants are independent of A and drop 
out of the ratio (|7.4|) . For Y C (A) and R 2 (lJ c ,A) see 
(|6T4|) and (|6T6|) . for # m (Y) see (j4~43l) . The anisotropy 
effect on the Binder cumulant t/(0; A) at T c will be de- 
scribed by the difference 



A£/ C (A) = C/(0;A) - [7(0; 1) 



where 



£7(0; A) = 1 - ^MYc(A)) [tf 2 (F c (A)] 2 



(7.7) 



(7.8) 
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VIII. QUANTITATIVE RESULTS AND 
PREDICTIONS 



For the application to three dimensions we shall em- 
ploy the same values as previously [H,[75|, A 3 — (47T)" 1 , 
v = 0.6335, u* = 0.0412, Q* = 0.945, B(u*) = 0.50. For 
reasons of consistency, a slightly different value will be 
used for a — 2 — Sv — 0.0995 in order to exactly satisfy 
the hyperscaling relation (|6.3[) . 



A. Universal bulk amplitude ratios 

Evaluating our analytic expressions for the bulk am- 
plitude ratios (|6.19|) and (|6.20|) in three dimensions we 
obtain for isotropic systems [compare (|3.12p ] 



f'l & = Qi 



-0.119, A-/A+ = 2.04 



.1) 



This can be compared with the series expansion re- 
sults for the three-dimensional Ising model by Liu and 
Fisher [75j who calculated the amplitude ratios (R^) 3 — 
0.0188 ± 0.0001 and A+/A~ = 0.523 ± 0.009. These cal- 
culations were carried out for several different cubic (sc, 
bcc, fee) lattice structures in order to test bulk univer- 
sality (see also [74|). The relation between Rf and Qi 



is (R^) 3 = — a(l — a)(2 — a)Q\. This yields the central 
values of the Ising model based on series expansions 



Q 



l\Is 



-0.1099, A-/A+ {I 



\Ising 



1.91 



Considering the fact that our present theory is an ef- 
fective finite-size theory that is not designed to produce 
highly accurate bulk predictions the results (|8.ip are in 
acceptable agreement with (|8.2[) . As seen from (I6.17P - 
(|6. 19|) , the bulk results for the free energy are sensitive to 
the choice of the geometrical factor in defining the renor- 
malized coupling, (|5.ip . The results (|8.ip demonstrate 
the appropriateness of the choice of Ad, (|5.2p . 
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FIG. 5: Finite-size amplitude Tc(l) , (fIU2j) . ([83}, of the free 
energy density of isotropic systems in a cubic geometry at T c 
for n = 1 in three dimensions (open circle), and at e = 1, 
(|8.7p (star), of the e expansion (|6.25[l , MC data for the three- 
dimensional Ising model on sc and bcc lattices [34l. |35|. 



for the finite-size amplitude .F c (l) of the free energy den- 
sity at T c [H 



^c(l) 



- 0.657 ±0.03 
- 0.643 ±0.04 



(sc lattice) 
(bcc lattice) 



(8.3) 



are consistent with the universality hypothesis. Subse- 
quently the more accurate MC result at T c was obtained 



^c(l) 



MC 



0.625 ± 0.005 



(sc lattice) (8.4) 



which is also consistent with 

In three dimensions the numerical values of the quan- 
tities Z c , Jo,/i, and I2 in our analytical result (|6.12p are 
l c = 2.042, J (/l,l) = 1.6430, h(%,l) = -4.1581, and 
7 2 (^,1) = -15.4032. This yields the theoretical predic- 
tion for the finite-size amplitude 



0.6315 , 



•5) 



in excellent agreement with the MC results (|8.3p and 
(Fig. 5) 



B. Finite-size free energy of isotropic systems 

1. Test of the d — 3 theory : amplitude at T c 

In order to test the reliability of our finite-size theory 
we first consider the isotropic case A = 1, = ^q + , 
where accurate MC data by Mon [34]; [H, [36| are avail- 
able. 

The first set of data was obtained for the three- 
dimensional Ising model with NN couplings on sc and 
bcc lattices. These systems have different values of T c 
and different correlation-length amplitudes £0+ but both 
belong to the subclass of (asymptotically) isotropic sys- 
tems with A = 1. Within the error bars, the MC results 



2. Epsilon expansion at T c 

For comparison we also evaluate the result of the e 
expansion (|6.25p . For isotropic systems (A = 1) the co- 
efficients in (|6.25p are well defined. The numerical values 
are 

/o(l) = - 0.3302, h(l) = -0.4218, (8.6) 

where 1 denotes the 4x4 unity matrix. For e = 1 the 
terms up to 0(e 1/2 ) of (|oT23)) yield 

J- C (l) e= i = - 0.7520, (8.7) 

which is in less good agreement with the MC results (|8.3p 
and (EH). (Fig. 5) 
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FIG. 6: Scaling functions f ex (x;l), (|6^2|) , (IrTTOl) , (Irl23|) 
for n = 1 (solid line) and ^(x;l), (^37)) . (^32]) . (|rl38)) 
for n = oo (dashed line) for the excess free energy density 
of isotropic systems in three dimensions as a function of the 
scaling variable x = t(L/£o+) 1/V • MC result (full circle) for 
the Ising model on a sc lattice [351 ] . 

3. Finite-size scaling function 

In three dimensions the numerical evaluation of the 
scaling function T ex (x, 1)^=3 as given in (|6.22[) . (|6.10p . 
and (|6.23[) yields the curve shown in Fig. 6 in the range 
—4.5 < x < 6 (solid curve). This range corresponds to 
the central finite-size regime (b) mentioned in Sect. IV 
A. A minimum with 

F ex {x min ;l) d =3 = -0.701, x min = -0.910 (8.8) 

exists slightly below T c . For the subclass of isotropic sys- 
tems within the (d = 3, n = 1) universality class both 
the position x m i n and the value T ex {x m i n \ l)d=3 are pre- 
dicted to be universal numbers. This can be tested by 
MC simulations for families of three-dimensional Ising 
models with A = 1, (e.g. on sc, fee, or bee lattices with 
isotropic interactions) in a cube with periodic boundary 
conditions. The nonuniversal differences of these models 
are predicted to be absorbable entirely in different values 
of £o + . In two dimensions such tests of universality for 
the critical Binder cumulant of isotropic systems at T c 
have been performed very recently by Selke [59| . 

Our analytical result for !F ex (x; l)d=3 is not applica- 
ble far outside the range of x shown in Fig. 6. In the 
limits x — > ±oo this result does not correctly describe 
the exponential decay to zero in the regimes (a) and (c) 
mentioned in Sect. IV A, as expected. 

For comparison we also present the exact result for the 
scaling function J-^(x; 1) in the large - n limit in three 
dimensions. For A = 1 and d = 3 the numerical solutions 
of feTiOl) and (|Q9"j) at T = T c are P c (l) = 1.946 and 

TZ (0; 1) = Too (0; 1) = -0.526. (8.9) 

In Fig. 6 the scaling function T^(x; 1), (|6.32p , is shown 
in three dimensions (dashed curve). Unlike the case n = 



I, T^(x; 1) does not have a minimum at finite x below 
T c but has a slow monotonic decrease towards a finite 
negative constant T^(~ oo; 1) = —3.18. Above T c it 
decays exponentially to zero (but not with the correct 
exponential form, see Sect. X.). 

C. Three-dimensional anisotropy 

In the following we present quantitative predictions for 
the nonuniversal effect of a non-cubic anisotropy on the 
finite-size scaling functions in three dimensions. We il- 
lustrate the three-dimensional anisotropy effects for the 
example of A and A given in (|2.14p and (|2.15p of Sect. 

II. In this example, A(w) depends only on the single 
anisotropy parameter w, (|2.16p , — \ < w < 1. The cru- 
cial anisotropy function is given by (|4.46j) which enters 
the functions J , h, and I 2 defined in (I4.44|) and (|4.45p . 

















' — \ 


n= 1 

n = infinity 




\ 






II \ \ 

'/ V 
J \\ 

1 V 






\ \ 

\\ 








-0.4 -0.2 


0.2 


0.4 


0.6 






FIG. 7: Differences AF c (A(w)), (|8A0)) . and AJ c " c , 00 (A(it;)), 
(fSTlj) . of the finite-size amplitudes (|67T2|) and (|6739)l of the 
free energy density of anisotropic systems with the reduced 
anisotropy matrix A(w), (|2.15[) . in a cubic geometry at T c for 
n = 1 (solid line) and n = oo (dashed line) in three dimensions 
as a function of the anisotropy parameter w, (|2.16|) . 

First we consider the anisotropy effect on the finite- 
size amplitude of the free energy density for n = 1 and 
n = oo at T = T c as described by the differences 

AJ r c (A(u))) = T{0; A(w)) - T(0; 1), (8.10) 

AT C!0o {A{w)) = f M (0;AW) - J-oo(0;l), (8.11) 

where J-(0;A) and ^^(OsA) are given by (|6.12p and 
()6.39p . This is shown in Fig. 7 in the range —0.45 < 
w < 0.80. The anisotropy effect is well pronounced 
for both positive and negative values of u>, with a non- 
negligible n dependence. For both n — 1 and n = oo 
two maxima of almost equal heights exist at w max — 
—0.333 and w max — 0.500, with A.T c ^ max = 0.0167 and 
A-Fcmax — 0.0165, respectively, for n = 1. [The slight 
difference of the two heights is presumably not a con- 
sequence of the approximations made for n = 1; such 
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FIG. 8: Difference AU c (A(w)), (fTf)) , of the Binder cumulant 
(17.8[) of anisotropic systems with the reduced anisotropy ma- 
trix A(w), (I2.15p . in a cubic geometry at T c for n = 1, d = 3 
as a function of the anisotropy parameter w, (I2.16|l . 



a difference exists also for the exact result for n = oo 
where A:F Ci00imaa , = 0.0202 and &F c ,oo,max = 0.0199 at 
Wmax = —0.333 and w max = 0.500, respectively.] At 
w = —0.45 and w = 0.80 we predict the larger negative 
values AT C = —0.069 and AT C = —0.079, respectively 
(not shown in Fig. 7). 

The corresponding anisotropy effect on the Binder cu- 
mulant for n = 1 at T c is shown in Fig. 8 as de- 
scribed by the difference AU c (A(w)), (J7T7J) [90j]. Figs. 
7 and 8 imply that the relative anisotropy effect on 
the free energy AJ r c (A(w))/^ r c (l) for n — 1 is consid- 
erably larger than that on the Binder cumulant. For 
the free energy it is predicted to be of 0(2.5%) at the 
maxima which may be detectable in future MC simula- 
tions of the three-dimensional Ising model. By contrast, 
the corresponding relative effect on the Binder cumulant 
AU c (A(w))/U(0; 1) for n = 1 is predicted to be only of 
O(0.6%) at the maxima. It becomes quite large, how- 
ever, in the regime w < —0.45, as shown in Fig. 1 of 
and in the regime w > 0.8. Previous MC simulations 
[42| of the anisotropic three-dimensional Ising model in 
the range —0.48 < w < are in disagreement with our 
results for the Binder cumulant ( see also [Hj]). 

In Fig. 9 we also show the predicted relative anisotropy 
effect on the susceptibilities \ + and X~ f° r n = 1 and on 
xto f° r n — 00 a t T c in the same range of w. While this 
effect is of 0(1%) near the maxima of the susceptibilities 
X + and xto the corresponding effect on the susceptibility 
X~ is only of O(0.1%). Previous MC simulations [I2,[43| 
on x~ did not resolve this small anisotropy effect. 

Our prediction of the anisotropy effect on the finite-size 
scaling function T ex (x\A(w)) for n = 1 near the mini- 
mum below T c is shown in Fig. 10 for several w. While 
the position of the minimum x m i n depends only weakly 
on the anisotropy the value T ex {x m in'-, 1) is significantly 
changed relative to the isotropic case (dotted curve in 
Fig. 10). This effect is well outside the error bars of the 
MC data by Mon for the isotropic case [34], [3f| and may 
be be detectable in future MC simulations. 



FIG. 9: Relative anisotropy effect on the susceptibilities x + 
and \~ f° r n = 1, (|7.2[> , and xto f° r n = °°! (|6,35[l . a t 
T c of anisotropic systems with the reduced anisotropy ma- 
trix A(w), (12.15[) . in three dimensions in a cubic geometry 
as a function of the anisotropy parameter w, (|2. 16|) . as de- 
scribed by Ax*(A(w)), (|7.4|) (solid and dot-dashed lines), 
and Axtc,c(A(w)), (fOB) (dashed line). 
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FIG. 10: Scaling function T ex (x; A(w)) , (^22j) . (pdO)) . (p^3|) 
of the excess free energy density of anisotropic systems with 
the reduced anisotropy matrix A(w), (I2.15|l . for n = 1 in 
a cubic geometry in three dimensions as a function of the 
scaling variable x — t(L' /£' 0+ ) 1 ^ 1 ' for several values of the 
anisotropy parameter w, (|2.16[) : w = 0.50, —0.33 (dot-dashed 
line), w = 0.80 (solid line), w = —0.45 (dashed line), w — 
(dotted line). MC result (full circle) for the three-dimensional 
Ising model on a sc lattice [35|. 



D. Two-dimensional anisotropy 

Highly precise numerical information on the nonuni- 
versal anisotropy effect on the critical Binder cumulant U 
of the two-dimensional Ising model has been provided re- 
cently by MC simulations of Seiko and Shchur [44( . They 
considered finite square lattices with isotropic ferromag- 
netic NN couplings K x = K y = K > and an anisotropic 
NNN coupling J only in the ±(1,1) directions but not in 
the ±(—1,1) directions (Fig. 11 (a) ). They found a non- 
monotonic dependence of U on the ratio J/K (as shown 
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in Fig. 4 of Ref. [4J). 

The anisotropy matrix of the corresponding two- 
dimensional if 4 lattice model is [H} 



A 2 =2a 2 K + J J 

\ J K + J 

with the reduced anisotropy matrix 

A 2 (s) = A 2 /(det A,) 1 / 2 = {l-s 2 )- 1 ' 2 

and with the single anisotropy parameter 
' =(l + K/J)- 1 




K 



J 



.12) 



(8.13) 



(8.14) 



By universality it is expected [13j that (in some range 
of A 2 near 1 and for sufficiently large L) the two- 
dimensional ip 4 model has the same anisotropy effects 
at T c as the two-dimensional Ising model if both models 
have the same reduced anisotropy matrix A 2 . Unfortu- 
nately, at the present time, it is not known how to per- 
form quantitative finite-size calculations for the ip model 
in two dimensions. 

It is possible, however, to incorporate a two- 
dimensional anisotropy of the type shown in Fig. 11 (a) 
in a three-dimensional ip 4 (or Ising) model on a simple- 
cubic lattice with isotropic NN couplings K x — K y — 
K > 0, with an anisotropic NNN coupling J\ = J ^ 
in the x — y planes, and with an additional NN coupling 
Kq > in the z direction (Fig. 11 (b)). The correspond- 
ing anisotropy matrix is 



2d 



7,2 



K + J J 
J K + J 
K 



.15) 



The eigenvalues and eigenvectors are Ai = 2a 2 (K + 2J), 
A 2 = 2h 2 K, A 3 = 25 2 X n , 




,(2) 




,e( 3 > = 






(8.16) 

The eigenvalues are positive in the range — | < J/K < 
oo, K > 0, K > 0. In the limit J/K — ► oo and J/K — > 
oo (or K — > + and K — > + at finite J > 0) the model 
represents a system of variables tpi on decoupled one- 
dimensional chains with ferromagnetic NN interactions. 

In its present form the matrix A3, (|8 . 15|) . contains both 
a two-dimensional anisotropy due to J and an additional 
anisotropy due to the coupling Kq. Although A3 con- 
tains A 2 as a decoupled 2x2 submatrix it is not ex- 
pected that, for general fixed Kq, this three-dimensional 
model exhibits the same type of anisotropy effect (as 
a function of the ratio J/K) as the two-dimensional 
model with the matrix (|8.12p . The reason is that it is 
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(b) 
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FIG. 11: Lattice points of the primitive cell of (a) a square 
lattice, (b) a simple-cubic lattice, with isotropic NN couplings 
K x = K y = K (solid lines) , with an anisotropic NNN coupling 
J in the x — y planes (dashed lines), and a NN coupling Ko in 
the z-direction (dotted lines). The corresponding anisotropy 
matrices A2, A2, and A3, A3, are given by (|8.12p . I|8.13[l . 
and (|8.15|) . (|8.17|) . respectively. 



not A 3 itself but rather the reduced anisotropy matrix 
A3 = A 3 /(det A3) 1 / 3 that governs the anisotropy effect 
according to the results of the preceding sections. This 
matrix is given by 



[*o(l - s 2 )} 



-1/3 



1 s 



\0 Kq 



Kq = 



Kq 



(K + J)=l< {1 - S) - 



.17) 



(8.18) 



We see that the J dependence of A3 differs qualita- 
tively from that of A 2 because of the additional s de- 
pendence of Kq at given Kq/ K > 0. What is needed is a 
kind of isotropic extension of the two-dimensional matrix 
(|8.13p to three dimensions parallel to the proposed four- 
dimensional extension (|6.29[) of the three-dimensional 
matrix A, (|6.28[) . This is achieved by the choice Kq = 1 
or Kq — K + J. Then the reduced anisotropy matrix 
becomes 



A 3 (s) = (1 



„2s-l/3 




(8.19) 



with the same anisotropy parameter s as in (|8. 13|) , (|8.14l) . 
[Naively one might have expected that the choice of the 
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additional NN coupling should have been K = K in- 
stead of K = K + J . But in this case the third 
diagonal element of A3 would become (A^) zz = (1 — 



2)-V3(i _ s )2/3 rat her than (1 - s 2 ) 



-13 



which would 



produce a qualitatively different anisotropy effect that is 
not an even function of s.] 

We have evaluated numerically the expressions (|7.8[) . 
(|4.43[) , (|6.14p in three dimensions for the Binder cumu- 
lant [/(0;A 3 (s)) at T c using the matrix A 3 (s), ([87T9]) . 
The result is shown in Fig. 12 in the range —0.8 < 
s < 0.8. The range < s < 0.8 corresponds to the 
range < J/K < 4.0 studied by Selke and Shchur [i3 |. 
The anisotropy effect shows up as a non-monotonic de- 
pendence on s. It is an even function of the anisotropy 
parameter s and exhibits two maxima of equal height 
at s max = ±0.461 corresponding to J/K = 0.855 and 
J/K = —0.316. This symmetry is a consequence of the 
symmetr y prop erty K 3 (y, A 3 (s)) = K 3 (y, A 3 (-s)) of the 
function (|4.46p . The symmetry is hidden if U is plotted 
as a function of J/K in which case the curve is asym- 
metrically distorted (see our Fig. 13 and Fig. 4 of [44|). 
Our theoretical value 17(0; 1) = J7(0;A 3 (0)) = 0.417 
for the isotropic three-dimensional </? 4 theory [66[ differs 
somewhat from the MC result [9l| 0.465 of the three- 
dimensional Ising model on a sc lattice and is, of course, 
far from the MC result 0.6107 [3!| of the two-dimensional 
Ising model on a square lattice. The magnitude of the 
anisotropy effect, however, turns out to be rather insensi- 
tive to the precise value of U (0; 1) of the isotropic system. 

To exhibit clearly the deviations from isotropy and for 
the purpose of a comparison with the MC data [44| for the 
anisotropic two-dimensional Ising model we have plot- 
ted in Fig. 13 our theoretical result for the difference 
AJ7 c (A 3 (s)), (|7.7p . as a function of J/K together with 
the corresponding difference of the MC data of Fig. 4 by 
Selke and Shchur [44J. The theoretical maximal value is 
AU c<ma x = 0.0010 at J/K = 0.855 and J/K = -0.316. 
The isotropic value, i.e. AU C — 0, is found at s = 
and at s = ±0.6169 corresponding to J/K = 1.611 and 
J/K = —0.3815. We see that for positive J/K there is re- 
markable agreement between the MC data in two dimen- 
sions and the anisotropic ip 4 theory in three dimensions, 
thus confirming our expectation regarding the similarity 
of the anisotropy effect in the two- and three-dimensional 
models. It should be noted, of course, that no exact 
agreement can be expected. Only the anisotropic two- 
dimensional ip 4 theory (with n = 1) is expected to yield 
exactly the same anisotropy effects (in the asymptotic 
critical region) as the two-dimensional Ising model. 

The non-monotonicity for small negative values of J/K 
and the maximum at J/K = —0.316 predicted by our 
theory was not detected in the preliminary MC simu- 
lations by Selke and Shchur [44j who found a mono- 
tonic decrease of U when taking a weak antiferromag- 
netic coupling J [92| . It would be interesting to perform 
more detailed MC simulations for the anisotropic two- 
dimensional Ising model in the regime of negative J/K. 

Very recently such MC simulations have been started 
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FIG. 12: Binder cumulant U(0; A 3 (s)) , |£Bj!, (TPS)) . (^Hl) , 
for anisotropic systems with the reduced anisotropy matrix 
A3(s), (18.19[) . in a cubic geometry at T c for n = 1 in three 
dimensions as a function of the anisotropy parameter s, (|8.14|l . 



by Selke [57f in order to test our prediction for the 
Binder cumulant in the regime of negative values of 
J/K. The positive value of his MC result ^ AU^ IC = 
0.00056 ± 0.00015 for J/K = -0.25 indeed confirms the 
predicted increase of AU C for small negative J/K accord- 
ing to Fig. 13. More quantitively, there is indeed reason- 
able agreement with our theoretical result AU C — 0.00073 
at J/K = —0.25 (as shown in Fig. 13) corresponding to 
s = —1/3. It remains to be seen whether the predicted 
symmetry with regard to s is also confirmed by MC sim- 
ulations. 



0.002 

-0.002 \- 
-0.004 
-0.006 
-0.008 
-0.01 
-0.012 



AU, 




— d=3 <]> theory 

* Kamieniarz and Bloete 1993 

a Selke and Shchur 2005 



2 

I/K 



FIG. 13: Difference A(7 c (A 3 (s)), fTf) . of the Binder cumu- 
lant f/(0; As(s)) shown in Fig. 12 but plotted as a function 
of J/K, together with the corresponding difference of the MC 
data of Fig. 4 of [3] for the two-dimensional Ising model. 



We also apply our general result (|6.12[) for the finite- 
size amplitude of the free energy density at T c to the case 
of the two-dimensional anisotropy determined by the ma- 
trix (|8.19[) . The anisotropy effect as described by the dif- 
ference AjF c (A 3 (s)) is shown in Fig. 14 for n = 1 (solid 
curve). The curve is an even function of s and has two 
maxima of equal height at s rnax = ±0.450 corresponding 
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to J/K = 0.818 and J/K = -0.310. The theoretical 
maximal value is AJ 7 c max — 0.0060 for n = 1. The cor- 
responding effect for n — oo (dashed curve) as computed 
from the exact result (|6.39[) is slightly more pronounced 
than for n = 1. Nevertheless the anisotropy effect for 
n = 1 may be detectable by MC calculations for both the 
three-dimensional and two-dimensional anisotropic Ising 
models. The two-dimensional model is of course a better 
candidate, as noted by Selke and Shchur because the 
value of T c is known analytically as a function of J/K . 
Although the solid curve in Fig. 14 is calculated on the 
basis of the ip 4 theory in three dimensions with the re- 
duced anisotropy matrix (|8.19[) we predict that this curve 
should be close to the difference AJ 7 c (A2(s)) of the free 
energy density of the two-dimensional Ising model with 
the reduced anisotropy matrix (|8.13p . It would be inter- 
esting to test this prediction by MC simulations. 

For comparison with Fig. 10 we have also computed 
the nonuniversal anisotropy effect on the finite-size scal- 
ing function of T ex for the reduced anisotropy matrix 
(|8.19[) near the minimum below T c as shown in Fig. 15 
for several values of s. Again this effect for s = ±0.80 
is well outside the error bars of the MC data by Mon 
34, 35] for the isotropic case and may be be detectable 
in future MC simulations. 



0.008 



0.006 



-0.002 - 



-0.004 - 





/ 

/ 

// 


\ 

\ 

, — V \ 


/ \ 
/ ^— » \ 




// 
// 


\ \ 
\ \ 


' / \\ 




// 


\\ 
\\ 


f\\ 




// 
« 


\\ 
\\ 








\\ 
\\ 


// \ 
// \ 








N\ 
\. 


y \ 


















n = 1 










n = infinity 










isotropic 


















-0.6 -0.4 



-0.2 
s = 



0.2 
( 1 +K/J)" 1 



0.4 0.6 



FIG. 14: Differences AT c (A 3 (s)), ([8^0), a nd A T CtOC ( A 3 (s)), 
(JSTTTJ), of the finite-size amplitudes (|Q2)) and (|Q9)) of the 
free energy density of anisotropic systems with the reduced 
anisotropy matrix A 3 (s), (18.19[) . in a cubic geometry at T c for 
n — 1 (solid line) and n = oo (dashed line) in three dimensions 
as a function of the anisotropy parameter s, (|8.14|l . 



E. Limit 



1 and Lifshitz point 



In Subsections C and D we have assumed the positivity 
of all eigenvalues X a ,a — 1,2,3. They vanish in the 
limits w — > 1 , w — ► — 1/2 and s — > ±1 in which cases 
our results for AU C and AT C are not applicable: they 
become singular as indicated by the curves in Figs. 7 - 
9 and 12 - 14. In the following we confine ourselves to a 
brief discussion of the limit s — > ±1 of the model shown 
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FIG. 15: Scaling function F ex (x; A 3 (s)) , (f630]> . 
of the excess free energy density of anisotropic systems with 
the reduced anisotropy matrix A 3 (s), (|8.19p . for n = 1 in 
a cubic geometry in three dimensions as a function of the 
scaling variable x — t(L'/^Q + ) 1// " for several values of the 
anisotropy parameter s, (|8.14|l : s = 0.45,-0.45 (solid line), 
s = 0.80,-0.80 (dot-dashed line), s = (dotted line). MC 
result (full circle) for the three-dimensional Ising model on a 
sc lattice 13511. 



in Fig. 11. (A similar discussion could be given for the 
model shown in Fig. 3 in the limits w — ► l,w — > —1/2.) 

(i) According to (|8 . 14[> . the limit s — > 1 can be per- 
formed as K — > 0+ at fixed J > and Kq > (keeping 
Ai > and A3 > positive while A2 — > 0+). In this limit 
our model is reduced to decoupled two-dimensional lat- 
tices which have ferromagnetic NN couplings J and Kq in 
the ±(1,1,0) directions and in the ±(0,0,1) directions, 
respectively. Such a model has a ferromagnetic critical 
point of the (d = 2, n = 1) universality class. Therefore 
it is expected at the outset that the results of our (p 4 
theory at fixed d — 3 must break down for s — > 1 . 

(ii) To discuss the limit s — > — 1 we first perform a 
rotation in wave- vector space, q = Uak, by means of the 
orthogonal matrix determined by the eigenvectors (|8.16p . 



U, 




(8.20) 



Correspondingly the inverse propagator ro±5i^(k) of the 
Hamiltonian (|2.7|) is transformed to r + 5K(q), where 
SK(q) = 5K(XJ^ 1 q) with the interaction part 



3 d 
6K(q) = ^ A «<7q + X! B &l3iS qaqpQfQS + 0(q 6 

a—l a,/3,7,<5 



.21) 



where A Q is given after (|8.15D . 

On the level of Landau theory a Lifshitz point exists 
at Ai = corresponding to s = — 1 or J/K = —1/2 
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with a wave-vector instability in the (1,1,0) direction. 
It is expected, however, that, due to fluctuations, the 
Lifshitz point occurs at a shifted value Ai = X%lp that 
depends on nonunivcrsal details of the model. [Similarly, 
roc depends on all details of the interaction, see (|4.32|) .] 
This corresponds to a shifted coupling ratio (J/K)lp = 
-1/2 + \i LP /(2a 2 ), presumably with Xilp < 0. To de- 
scribe the critical behavior near the Lifshitz point would 
require to introduce a renormalized shifted eigenvalue ac- 
cording to Xr = Z\(Xi — Xilp) [§3|- It would be interest- 
ing to locate this Lifshitz point by MC simulations and to 
detect the nonuniversal change of finite-size effects [crit- 
ical Binder cumulant and free energy at T C (J/K)] upon 
approaching this point along the "A-line" T = T C ( J/K) 
as J/K — ► (J/K)lp. This change can be compared with 
our predictions shown in the curves of Figs. 12 - 14 for 
negative J/K and negative s. 



IX. HYPOTHESIS OF RESTRICTED 
UNIVERSALITY 



The results for the finite-size scaling functions (| 6 . L [) 
- (|6l6]l . (pL22]) . (jQ2"l) . ([731) . and (HH) depend on the 
nonuniversal anisotropy matrix A but are independent 
of the bare coupling ito, of the lattice spacing a, of the 
cutoff of Lp A field theory, and of the fourth-order moments 
Ba/3-yS etc.. We anticipate that the finite-size scaling 
functions would also remain independent of higher-order 
couplings, such as those of ip 6 terms and of higher-order 
gradient terms etc., if they were included in the Hamil- 
tonian. A special matrix A with given matrix elements 
A a/ 3 can be obtained from various different lattice struc- 
tures with a large variety of different couplings, both in 
0(n) symmetric ip 4 lattice models and in 0(n) symmet- 
ric fixed-length spin models. We expect that J-(i, A) and 
U(x 1 A) are the same for all those systems whose geom- 
etry and boundary conditions are the same and whose 
reduced anisotropy matrix A is the same. We consider 
this feature as a kind of restricted universality within a 
(d, n) universality class. A nontrivial aspect of this fea- 
ture is that it is governed by the bare anisotropy matrix 
A containing the unrenormalized microscopic couplings. 
Our approximate results do not yet provide a rigorous 
proof for the validity of this hypothesis. It would be in- 
teresting to test this hypothesis by MC simulations for 
microscopic spin models with such anisotropy matrices. 

For concreteness consider the three-dimensional 
anisotropic model (i) of Sect. II A with the three cou- 
plings K, J, and K as described by the matrix 



D J + K J + K 
A = 'l'<r | J + K D J + K 
J+K J+K D 



(9.1) 



with D = K + 2 J + K and with the reduced matrix A, 
(I2.15p . For a given fixed value of the anisotropy param- 
eter w, (|2.16p . a family of anisotropic spin models with 



different couplings K, J, and K are predicted to have the 
same finite-size scaling functions if the third-NN coupling 
K is chosen as 



K 



1 



1 



[wK- (1 - 2w)J] 



(9.2) 



in the range where A Q > 0. Equation (|9.2[) represents a 
surface in the space of the three couplings K, J, and K. 
At K = 0, this surface becomes a "A - line" 



J 



l-2w 



-K 



(9.3) 



along which, at a given fixed value ofu>, —1/2 < w < 
1/2, all finite-size scaling functions of models with the 
anisotropy matrix (|9.1j) are predicted to remain un- 
changed when changing K and J simultaneously accord- 
ing to ([92]) (in the range K + J > 0, K + AJ > 0). 

A non-trivial test of our hypothesis applied to two di- 
mensions can be performed for the following example. 
Consider a triangular-lattice model with a shape of a 
rhombus and with three NN couplings K\, K2, K 3 , and 
a NNN coupling J only in the ±(3/2, V3/2) directions 
(Fig. 16). The anisotropy matrix of this system with 
lattice constant a = 1 is 



AK-y + K 2 + K 3 + 9J \/3 (K 2 - K 3 + 3 J) 
V3 (K 2 -K 3 + 3 J) 3(K 2 + K 3 + J) 



(9 



In the absence of the NNN coupling J, isotropy is possible 
only for the symmetric case K\ = K 2 = K 3 . In this case 
the critical Binder cumulant for the (n = l,d = 2) uni- 
versality class for periodic boundary conditions is known 
to very high accuracy [Hj]: U = 0.6118277 ± 0.0000001. 
Apart from this case, the system can become isotropic 
even for K 3 ^ K = K\ = K 2 if the NNN coupling J is 
chosen as 



according to 



J = \ (K 3 - K) 



in which case 



(9.5) 



2K + AK 3 



2K + AK 3 ) 



(9.6) 



and A = 1. Then, on the basis of our hypothesis 
of restricted universality, the critical Binder cumulant 
for K 3 / if, J ^ is predicted to have exactly the 
same value as found by Kamieniarz and Blote [391 ] for 
J = 0, K 3 — K . A corresponding prediction should hold 
also for other boundary conditions (e.g., free boundary 
conditions). 
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FIG. 16: Lattice points of a triangular lattice with the shape 
of a rhombus. The solid lines indicate the anisotropic NN cou- 
plings Ki, K2, Kz, the dashed line indicates the anisotropic 
NNN coupling J. The anisotropy matrix A is given by 



X. 



FINITE-SIZE EFFECTS OUTSIDE THE 
CENTRAL FINITE-SIZE REGIME 



The result F ex {x; A), (pT22)) . (|6J0]) . needs to be com- 
plemented outside the central finite-size regime where 
this result does not have the correct exponential structure 
for large \x\. This is the regime below the dashed lines 
in Fig. 1. It turns out that it is necessary to further dis- 
tinguish between a scaling and a nonscaling regime (the 
latter is the shaded region in Fig. 1). Both regimes be- 
long to the asymptotic critical region. In these regimes 
ordinary perturbation theory is appropriate. Here we 
perform the corresponding analysis above and below T c 
at the one-loop level. In order to distinguish the pertur- 
bation results of this section from those of the preceding 
sections we explicitly include the indices + and — in the 
notation of T ex '+ , /'+, T ex -~ , and /'~. 



A. Scaling regime 

The starting point of ordinary perturbation theory for 
the bare free energy density (|2.25|) for n = 1 at h' = is 
in one-loop order [i.e., up to 0(1)] (see App. B) 

/'+ = ti(r ,L',K itj ,v')+O(u' ) , (10.1) 



= 2 r " M ™/ + u Xf + /o(- 2r °' L', K itj ,v') + O(u' ) 

(10 

above and below T c , respectively, where is given by 
flUBt and 

t ' ts ../^ M**) 



f- = -^-r[-2(ro-r 0c )] 2 + |(ro-r 0c ) j -J— 
64«o 2 J SK'(k' 



+ fo(-2(ro-r 0c ),L',K t j,v'). 



(k') 
(10.5) 



We define the finite-size parts Sf'^ of / /=t in the same 
way as 5f in (|4.33p and (|4.34p . Calculating the sum in 
the continuum limit v' — ► at fixed \ro — rod (see 
App. B and C) one obtains for 2 < d < 4 

8f' + (ro - r 0c , u' , L>, A) = ~ (r - r 0c ) d / 2 
+ Go((r - r 0c )L' 2 : A) + O(u' ) , (10.6) 



5f (r -r Qc ,u' ,L',A) 



64u 



7 [-2(r - r 0c )Y 



-^[-2(r -roc)] d/2 + Go(-2(r - r Qc )L> 2 ; A) 



+ O{u' ) 



(10.7) 



where G Q is given by ([B1H]) . Eqs. pU7B|) and ([HTTP cor- 
respond to (|4.37[) . The renormalized counterparts of 
<5/ /:t: are defined in the same way as f' R in (|5.3[) . (|5. 17|) . 
The explicit form of the functions f' R depends on the 
choice of the flow parameter. For the application to the 
regime \x\ ^> 1, we make the bulk choice = r(l+) for 
T>T C and n 2 P_ = -2r(Z_) for T < T c , with fi' 1 = £ 0+ , 
(|5.16p (the same choice will be made for the calculations 
in subsection B below). Then the functions f' R are given 
by 

f'+(r(l + ),u'(l + ),l + n,L',A) = - A d (l +f i) d /(id) 
+ ^TM l l^ L ' 2 ^ A ) + W+)) > ( 10 - 8 ) 



/^(r(?_), U 'G_),L M ,L',A) = 

+ 0(u'(i_)). (10.9) 



For l± — ► 0, this leads to the finite-size scaling function 
2)of the excess free energy density in one-loop order in the 
limit of zero lattice spacing 



f (r ,L ,K itj ,v ) = 



^]>>{[r„ + ^'(k')]K) 2/d }. (10.3) 



2L' d 



Because of the k' = term, the sum exists only for tq > 
0. Rewriting these expressions in terms of ro — ro c with 
roc given by (|4.32p yields up to O(l) 



where — £,' 0+ t~ 1 ' and 



0{u*) ,(10.10) 



+ 0{u*) (10.11) 



are the bulk second-moment correlation lengths above 
and below T c , respectively (for £q + see (|5.16p ). Here we 
have confined ourselves to the simplest form of pertur- 
bation theory for the regime \x\ ^> l,x = t{L' /£ QJ V~I V . 



f' + = f Q (r -r Qc ,L',K ZJ ,v'), 



(10.4) As a shortcoming of this approach, T^J^ diverges for 
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x — > at fixed finite L' which originates from the k' = 
term of (|10.3p . [This divergence could formally be sup- 
pressed by an L 1 dependent choice of l± but this would 
not avoid a structurally incorrect nonanalytic t depen- 
dence at t = for finite L'.] 

^l^laopix; A) serves the purpose of complementing 
F ex (x;A), (|fT22"l) . (j6T0)) in the large x regime. This 
is illustrated by the thin solid line in Fig. 17 for the 
example of three-dimensional isotropic systems (and for 
systems with cubic symmetry) with A = 1, V = L, 
£o+ = £o+ , Co- = Co- • The curves match reason- 
ably well above T c . [No perfect matching can be ex- 
pected because of the missing 0(u*) terms in (|10.10p 
and because F ex (x;l) is not applicable to the region 
i > 1 where it has an algebraic approach to a finite 
limit ^(oo; 1) = -2u* = -0.082 for x -> oo.] 

By contrast, ^il'i^ op (x; A) and T ex {x;A) do not 
match well below T c for two reasons: (i) The two- loop 
terms of 0{u*) in (|10.10[) are non-negligible, (ii) our 
approximate result !F ex (x;A) as represented by (|6.22|) . 
(I6.10p . (|6.23p is not applicable to the region x < —5. 
[In this region this result for T ex (x\ 1) has an unphysical 
maximum T ex (x ma x\ 1) = —0.303 at x max = —5.61 and 
has an algebraic approach to a finite limit T ex {— oo; 1) = 
—0.49 for x — > — oo.] Thus substantial further work is 
needed for a satisfactory description of the region well 
below T c . 

Nevertheless, Fx^^op nas ^ ne advantage of displaying 
the expected exponential large \x\ behavior. The leading 
large x behavior of Go for the isotropic case is (see App. 
B) 



/ x v \ (rf-l)/2 

-2d[ — ) exp(-i l/ ) + 0(exp(-2i l/ )). 

(10.12) 



For d = 3, Eqs. (fT0~T2l) and (fTOTOl) yield for \x\ v > 1 
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FIG. 17: Scaling functions T ex {x;l), (TQ2j>, (IBTTOj) . (TQ3|) 
for d = 3 (thick solid line), J?[l'£ op (x; 1), (|10.10|) for d = 3 
(thin solid lines), and (I10.13|l (dashed lines) for the excess 
free energy density of isotropic systems as a function of the 
scaling variable x = t(L/^o+) 1 ^ . MC result (full circle) for 
the Ising model on a sc lattice [35|. No scaling function ex- 
ists in the large - \x\ regions above and below T c which are 
sensitive to all nonuniversal details of the model according to 

^^(L/fcfcsija/et), BEEP , with 



B. Non-scaling regime 

1. Anisotropic ip 4 lattice model with finite lattice constant 

So far we have taken the continuum limit which is 
well justified in the range shown in Figs. 10, 15, and 
17 provided that £'±/a ^> 1 is sufficiently large. In 
earlier work [9| it was pointed out for the example of 
the susceptibility that the finite lattice constant a be- 
comes non- negligible in the limit of large L/a at fixed 
T =/= T c in the regime where the finite-size scaling func- 
tion has an exponential form. Here we further discuss 
this issue in the context of the excess free energy of the 
model (|2.1[) with V = L d and cubic anisotropy, i.e., on 
a simple-cubic lattice with lattice constant a and only 
NN couplings K. In this case we have A = 2a 2 Kl, 
A = 1,L = (2a 2 K) 1 / 2 L / , a = (2a 2 K ) 1/2 a' , and there ex- 
ist well defined bulk second-moment correlation lengths 



as shown by the dashed lines in Fig. 17 (with = 
=2"|£n. 

In case of noncubic anisotropy, all curves in Fig. 17 
including the thin solid lines and dashed lines are, of 
course, affected by the anisotropy matrix A ^ 1 in a way 
similar to that shown in Figs. 10 and 15. It would be 
straightforward to illustrate this effect by complement- 
ing Fig. 17 accordingly by means of curves representing 
T e x {x- A 3 (s)) and ^'^ op (x; A s (s)) , with A 3 (s) given 
by (|8.19[) , for several examples of s. In this case the scal- 
ing argument (horizontal axis of Fig. 17) needs to be 
replaced by t(L' /£ 0+ ) 1/u . 



(10.13) f± = (25 2 K) 1 / 2 ^' ± above and below T c (for n = 1). As 
shown in App. B, the excess free energy density in one- 
loop order attains the following form in the limit of large 



L/a = L'/a' at fixed arbitrary a/£± = a'/C± > 



L-'F**(L/Z±;l;a/t ± ) (10.14) 



F ex ' ± (L/Z ± ;l;a/Z ± ) = -d 



1 



2£d 



d-l 



L 



2<± 



(d-l)/2 



ex Pi -- 



2£± ..fa 
— — arsmh — — 
a \2£± 



(10.15) 
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This result applies to the shaded region of Fig. 1. 
The exponential part of (|10.15[) can be rewritten as 
exp(— L/£ e ±) with the exponential correlation lengths 



e 



e± 



(10.16) 



above and below T c , respectively. As a nontrivial rela- 
tion between bulk properties and finite-size effects [l(|, 
the lengths £ e ± describe the exponential part of the bulk 
order-parameter correlation function [76j in the large- 
distance limit in the direction of one of the cubic axes at 
arbitrary fixed T ^ T c above and below T c (for n = 1), 
respectively. This relation is exact in the large-n limit 
above T c 0. 

It has been shown Q that, because of the exponen- 
tial structure of the finite-size part of the susceptibility, 
the a dependence of £ e ± cannot be neglected even for 
small a/£± < 1 if L/£± >j241n2](£ ± /a) 2 is sufficiently 
large (see Fig. 1 of Ref. The same argument now 

applies to the a dependence of the exponential part of 
re*,±(L/Z±;l;a/t,±), ([TOTTS) . This implies that finite- 
size scaling and universality are violated in the large - \x\ 
tails oi J- ex,± at any a/£± > even arbitrarily close to T c 
because ultimately, for \x\ — > oo (i.e., for large L at fixed 
\t\ 0), the tails of ^ rea; ' ± become explicitly dependent 
on a. (Below we shall show that the tails depend also on 
the bare four-point coupling uq.) Thus no finite-scaling 
form (|1.2p . (|1.3[) . or (|1.4|) with a single scaling argument 
oc tL 1 l u and with a single nonuniversal amplitude C\ 
can be defined in this large - |i| region (94J. Higher- loop 
contributions cannot remedy this violation. It is obvious 
that an even larger variety of different nonscaling effects 
exist in the exponential finite-size region of systems with 
non-cubic anisotropics (A ^ 1). 



2. Isotropic (p 4 field theory 

The diversity of nonuniversal non-scaling effects in the 
region L/£± ^> 1 discussed above exists not only in 
anisotropic lattice models but also in fully isotropic sys- 
tems. We demonstrate this point for the isotropic (p 
field theory based on the standard Hamiltonian 



H 



field 



J d d x[^+ l -{W V f + u^} (10.17) 

v 



in a cube with V — L d and periodic b.c. and with some 
cutoff A in k space. Keeping the cutoff finite may be 
a valuable tool for testing universality as has been con- 
vincingly demonstrated by Nicoll and Albright [95[ in the 
context of bulk universality [8J . In a similar spirit this was 
done in [t| [ll[ with regard to finite-size universality in 
the large-n limit at and above T c . We shall show that 
jrex,± depend on the bare coupling uq and on the cutoff 
procedure for large L above and below T c of the n = 1 
universality classes. For the case of a sharp cutoff we 



shall also correct a previous misinterpretation [llj, |96j of 
the singular part of the excess free energy density at T c . 

Since f ex '^ has a finite limit for A — * oo we first we 
calculate f ex <^ at infinite cutoff A = oo within the mini- 
mal renormalization scheme at fixed dimension 2 < d < 4 
[62l . I&3L |82| . In one- loop order we obtain 



where for large L/£± 



L 



(d-l)/2 



(10.18) 



27r£± 



exp<i-— |(10.19) 



with the bulk second-moment correlation lengths 

C±(t;u) = ^(u)^ {I + C ± (t,u)} . (10.20) 

(There is no difference between exponential and second- 
moment correlation lengths at infinite cutoff at the one- 
loop level.) The function C±(t,u) represents the Wegner 
series 



C±(t,u)= 5>i m) (u) \t\ Ar ' 



(10.21) 



with the universal Wegner exponent A = lov, lo = 
d(3 u (u, e)/du\ u=U i,, and the Wegner amplitudes a^\u) 
depending in the nonuniversal renormalized coupling u. 
The latter is defined by 



u = A d Z u (u,e) 1 Z v (u,e) 2 u i 



0+ 



(10.22) 



(with the choice [i = £nl) where Z u {u, e) and Z v {u, e) are 
the standard Z factors [86| . Equation (|10.22p determines 
u as an implicit function of wo£o+- Although C± is an 
negligible additive correction in (110.20(1 for sufficiently 
small \t\ this is not the case in the exponential part of 
(110. 19)1 which, for small C±, can be rewritten as 



exp 



L 



A(L, t, u) exp 



(10.23) 



A(L, t, u) = exp < C±(t, u) 



L 



&±|t|" 



0(Ci) K10.24) 



with the nonuniversal nonscaling prefactor A(L, t, u) that 
cannot simply be replaced by 1 for small \t\. Even for 
arbitrarily small |t| ^ the prefactor becomes nonneg- 
ligible if L is sufficiently large, L ^> \C± | _1 ^o± \t\~ v ■ 
Thus the tails of the large L dependence of be- 
come nonuniversal and have a nonscaling L dependence 
through the prefactor A(L, t, u). This applies to the 
shaded area of the asymptotic critical region above and 
below T c shown in Fig. 1 . The same argument applies to 
the preceding subsection: it is necessary to keep the com- 
plete non-asymptotic form of the second-moment bulk 
correlation lengths at finite a 



U(t;u d e ) = £o±\t\- v {1 + C±{t,u a £ )} 



(10.25) 
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in (|10.15[) and (|10.16[) and to include all correction terms 
in C±{t,u a £ ). 94} 

The reasoning described above must also be extended 
to the case when a smooth cutoff A in k space is taken 
into account. This can be done by including an isotropic 
(Pauli-Villars type) term \(V 2 ip) 2 /k 2 in the Hamilto- 
nian (| 1 . 1 T|) [§, [9?j|. In this case the structure of F ex -± 
still remains exponential oc exp[— L/£ e ±(A)] for large 
L/£± but the exponential correlation lengths 

£ e± (A) = C±[l - \ A- 2 C± 2 + 0(A- 4 C ± 4 )] (10.26) 

become cutoff dependent. This causes a cutoff dependent 
prefactor A(L, t, u, A) in (|10.23p . 

As pointed out in fioj there exists a close relation be- 
tween the L dependence of finite-size effects and the x 
dependence of the bulk order-parameter correlation func- 
tion Gb discussed in Sect. III. In retrospect, the argu- 
ments presented above apply also to the exponential part 
of Gb oc |x|~ d+2 exp(— |x|/£ e ±) even if it is isotropic be- 
cause here the same correlation lengths £ e ± appear as in 
the large L decay of the finite-size quantities. No scal- 
ing functions <I>±, (|3.5|) , (|3.6p . (|3.19|) can be defined in the 
exponential large-distance regime |x|/£± 3> 1 (shaded re- 
gion in Fig. 2). Thus the exponential tails of G&(x; t) of 
the ip 4 theory have a nonscaling form that depends on uo 
and the (smooth) cutoff even for arbitrarily small t ^ 0, 
h = and t = 0, h 0. In addition, for anisotropic 
systems, it depends on the anisotropy matrix A and the 
higher order tensors B etc.. 

Although the nonscaling effect on the relative quantity 
j?ex,±jj7± becomes arbitrarily large for sufficiently large 
_L/£± this happens in a region where the magnitude of 
jrex,± j^ se if j s exponentially small. Thus, from a purely 
quantitative point of view, this is only a very small effect 
for systems with short-range interactions and periodic 
boundary conditions. 

This is in contrast to the corresponding non-scaling 
finite-size effects in the presence of (effective) long-range 
correlations caused by a sharp momentum cutoff —A < 
k a < A used in [{J [yj. Such a cutoff has often been 
used in the formulation and application of the RG the- 
ory based on the ip 4 Hamiltonian (| 10. 1 T[) (see, e.g., 
(9a , [H, [9§] ). As far as thermodynamic bulk properties 
are concerned this is well justified as the sharp cutoff docs 
not affect the critical exponents and the thermodynamic 
bulk scaling functions. Thus the ip 4 model (|10. 17|) with a 
sharp cutoff is a legitimate model of statistical mechan- 
ics that belongs to the same (d, n) universality class as 
systems with short-range interactions or with subleading 
long-range interactions. This implies the validity of ther- 
modynamic two-scale factor universality in the presence 
of a sharp cutoff. Chen and the present author 0, Hlj 
have raised the question whether this remains true also 
for confined systems. It was found, for the susceptibility 
and for the excess free energy in the large-n limit above 
T c , that a sharp cutoff is not compatible with an exponen- 
tial size dependence and violates finite-size scaling in the 



large - L regime above T r . This behavior was traced back 
to the well known [ll|) HH artifact that the sharp cutoff 
in k space causes long-range correlations in real space 
as can be demonstrated in the bulk order-parameter cor- 
relation function G{,(x;t; A) [ill \7l\ whose algebraically 
decaying non-scaling part dominates the exponentially 
decaying scaling part. By means of a RG one- loop cal- 
culation for n — 1 we find that this property holds both 
above and below T r for sufficiently large L. 

In contrast to [ll| . however, we do not obtain a vi- 
olation of finite-size scaling in the central finite-size re- 
gion including T — T c . Our present analysis is based on 
an appropriate decomposition of the excess free energy 
into singular and nonsingular parts in the sense of (|2.5|) 
whereas in [Tlj no L dependent nonsingular part was de- 
fined. We find that, in the presence of a sharp cutoff 
A and for large LA, Eq. (I10.18|) with (|10.19[) is to be 
replaced by 

/f =r A x f(t,L) + fZ ls (L), (10.27) 

with the singular part 

f?f(t, L) = L-^-^^A- 1 ) + /££(t, L), (10.28) 

poo 

* d (z) = / dy [e- (1+ ^y -e-y]E d (y), (10.29) 
Jo 



E d (y) = 




and the L dependent nonsingular part 

POO 

fA% s (L) =L- 2 A d - 2 dye-yE d (y). (10.31) 
Jo 

Although our one-loop result (jl0.27j) - (|10.3ip for the to- 
tal excess free energy den sity f^' (t, L) is equivalent to 
equations (8) and (16) of [1JJ, there is a crucial difference 
with regard to singular part. In contrast to the non- 
vanishing function &d,d'(z) of [TTj] , our function $d(z) 
vanishes at criticality, <&d(0) = 0. The temperature in- 
dependent part (|10.3ip cx L~ 2 should not be attributed 
to the singular part as was done in [Til ] . Our definition 
of the nonsingular part f^_ ns {L) oc L~ 2 is parallel to the 
standard analysis of bulk systems with a specific heat 
C± = A±\t\~ a + Cb with a negative critical exponent 
a whose finite value Cb at the finite cusp must not be 
included in the singular scaling part ~ but rather 

in the nonsingular " background" contribution of the spe- 
cific heat. The nonuniversal power-law term oc L~ 2 A d ~ 2 
in p0.28p dominates in the shaded region of Fig. 1 com- 
pared to the scaling part oc L ~ d but vanishes at 
T = T c and is subleading in the central finite-size regime. 
Thus the leading finite-size contributions in the <p 4 model 
with a sharp cutoff are in agreement with universal finite- 
size scaling in the central finite-size regime if the singu- 
lar part of the free energy is identified correctly. Conse- 
quently, the leading singular part of the Casimir force (in 
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film geometry) at bulk T c [llj remains universal within 
the subclass of isotropic systems even in the presence of 
a sharp cutoff but an additional regular part cx L~ 2 ex- 
ists that is nonuniversal and is dominant compared to 
the singular part cx L~ d . This unusual behavior is due 
to th e long-range correlations caused by the sharp cutoff 
[lOOl ] , as noted already in [ll| , which is of course a math- 
ematical artifact and not generic for systems with purely 
short-range interactions. As pointed out by Dantchev 
et al. [96j |. the sharp cutoff implies an unphysical dis- 
continuity of the slope of the interaction 8K{k) = k 2 at 
the boundary of the Brillouin zone which is the mathe- 
matical origin of the L~ 2 terms. For the reasons given 
above, however, we disagree with the opinion expressed 
in (96| that the concept of finite-size scaling as devel- 
oped for systems with short range interactions does not 
apply to the y? 4 model (|10.17p with a sharp cutoff. The 
authors of [96[ did not perform an analysis based on a de- 
composition of the type (|2.5| and (|10.27|) . Our analysis 
shows that, in spite of the mathematical artifact of K(k) 
at the Brillouin-zone boundary, the concept of finite-size 
scaling is well applicable to the central finite-size regime 
including T = T c and that a violation of finite-size scal- 
ing occurs only in the large-!/ regime at T ^ T c (shaded 
region in Fig.l), as in the other cases discussed above in 
the presence of a lattice cutoff or a smooth cutoff. 

Finally we discuss the case of an additional subleading 
long-range interaction of the van der Waals type as de- 
fined in (|2.7p and (|2.9|) . It was pointed out by Dantchev 
and Rudnick [l5| that it affects the finite-size susceptibil- 
ity in the regime L/£+ 1, similar to the effect caused 
by a sharp cutoff 0] . The effect of this interaction on the 
excess free energy f£ x and on the critical Casimir force in 
the case of film geometry was first studied by Chen and 
the present author [IJ, [lj| . The asymptotic structure for 
^> I in one-loop order above T c at h = is 



T ex (L/t; + )+bL 2 -°iS>(L/Z + ) (10.32) 



which is similar to (|3.17p . We have verified that, for 
n = 1, the same structure is valid also for cubic geometry 
with periodic b.c. above and below T c where the func- 
tion ^cube has an algebraic large-L behavior ~ (L/£±) , 
The latter dominates the exponentially decaying scaling 
part J? ex >± ~ exp(— L/£ e ±) m the shaded region of Fig. 
1. This implies that, in this region, two nonuniversal 
length scales &VO-2) an d £ ± a t h = govern the leading 
singular part of the excess free energy density 



L~ 



d 




cr-2 






. L - 




. L . 



(10.33) 



even arbitrarily close to criticality. In addition, there 
is the nonuniversal uq - dependent exponential tail of 
j7ez,± j n (|1Q.33|) . both the amplitude ~ b and the power 
—d — a of the L dependence are nonuniversal. Thus the 
universal scaling form (|1.2p . with only one length scale 
Cy V at h = 0, is not valid in the entire range of its scaling 



arguments for isotropic systems with van der Waals type 
interactions although such systems are members of the 
same (d, n = 1) universality class as, e.g., Ising models 
with short-range interactions. The structure of (|10.32[) 
and (|10.33p and the corresponding structure for the crit- 
ical Casimir force in film geometry has been confirmed 
in SlUIii. 
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Appendix A : Universal bulk amplitude relations 

In this Appendix we present explicit expressions for 
the universal constants Qi, Q3, and Pi, (|3.9p - (|3.15|) . in 
terms of universal scaling functions. Near T c the sum 
rule (j3~2oT) yields 

x' b (t,hf) = -A! 1 \t\*'&w ± (Ai l ti\t\-e s )/dh' 2 

= D[t' ± (t,ti) 2 -*>$ ± (D' 2 h'\t\-e s ) (A.l) 

with the universal function 

00 

$±(y) = 27: d / 2 T(d/2)- 1 J rfss 1 -' ) $ ± ( S , y) . (A.2) 



At t > 0,ti = 0, JA3|) yields x fc M) = F' + \t\-~< with 
T' + = -A[ A' 2 2 W 2 = D[ (^ + ) 2 -" $+(0) (A.3) 

where W 2 = lim^ d 2 W+(y)/dy 2 . At t = , h! f we 
have from 02) &(0, b!) = C h = CW\~ v/m with 

t' c = e a+ (D' 2 )-»/^X , (A.4) 

thus (TP) yields ?&(0,h') = VAh'^'^ with 



-A\ A' 2 1+1/5 W 



X 



2-r, 



= D[ [£o+ D 2 
where $(00) = $±(00) and 

W = lim {\yV /m d 2 W±{y)/dy 2 } 



$(00) (A.5) 



X = lim \y 

y-*oo 



{\y\ v/m X ± (y)} 

(|A~3|l and (|A~5|l yield 

r;/r c = a'^- 1 ' 5 w 2 w- 1 



(A.6) 
(A.7) 



D't 1/& $+(0) \X 2 ^ *(oo)l , (A.8) 



3G 



thus we obtain the universal ratio 



A' 2 

Pi = — 
D' 



W $+(0) 



W 2 X 2 -n $(oo) 



5/(5-1) 



(A.9) 



Eqs. (| A.4|) and (|A.8|) yield a universal ratio Q 2 different 
from P 2 , 



Q2 = 



2-v 



<MQ) 

$(00) 



(A.10) 



where we have used 1 — 1/5 = j/(0S) and (2 — r\)v = 7. 

Following Privman and Fisher [5| we assume that the 
unsubtracted bulk correlation function G^x^ — x^ ; t, h') = 
liniv^oo ( t / 5 '( x i) ( ^'( x j '))' nas the asymptotic scaling form 

G' b (x';t,ti) = D[\ X >\- d+2 -" Z ± (|x'|/£' ,D' 2 h'\t\-^) 

(A.ll) 

with the same constants D[ and D' 2 as in Eq. (|3 - 5[) and 
with a universal scaling function Z±(x, y). From (jA.lljl . 
(|2.3ip , and (|3.3|) we obtain the square of the bulk order 
parameter below T c 

[m' h (t)] 2 = lim lim G' h hc';t,h') 

|x'|->oo 

= D'i {^_)- d+2 -'i\t\^ d - 2+ ^ Z 
= lim [dm, h')/dhf = K A' 2 f\t\^ Wi (A.12) 

with Z = lim^o lim^oo {x~ d+2 ~ ri z_(x, y)} and Wi = 
lirrij,^o dW-(y)/dy. In order to derive Qi and -P3 we use 
(|A.3|) and (jA.12|) and obtain 

D[ = -A[ A' 2 2 (e 0+ r 2+ " W 2 /$+(0) 

= K A' 2 ) 2 (^) d - 2+ " Wi/Z . (A.13) 
Together with (|3.16p this yields the universal quantities 

Qi = A (e 0+ ) d = -z w 2 [x_(i))]- d+2 -^ 1 \w x $ + (o)" 

(A.14) 



and 



P 3 = D^A'-^'^iA^)- 2 = - Q- 1 '^' W 2 /$+(0) . 

(A.15) 

Finally we consider the universal ratio (|3.14[) . The am- 
plitude D'^ is given by D 1 = D[ $±(0,0) 6 with the 
universal constant [13, llOll 



)-7/(^) 



6 - DL 



$±(0,0) _ (47r) d / 2 r(^) 



$±(0,0) 2 d - 2 + r iT( d ~ 2 2 + ' n ) 



.(A.16) 



Together with (|A.3j) this yields a universal ratio Q3 dif- 
ferent from P3 , 

Q 3 = DUti+) 2 -yr> + = $ ± (o,o)c/i+(o). (A.17) 

The universal constant Q3 in (|3 . 1 5[) is 

Q 3 = $±(0,0)/$+(0) . (A.18) 



Appendix B : Gaussian model with lattice 
anisotropy 



In order to derive the Gaussian part of (j4.26|) and the 
results of Sect. X we consider the Hamiltonian (|2.1|) and 
(12. 7[) for ro — a$t > 0, uq = and h = with N scalar 
variables ipj on a simple-cubic lattice with lattice con- 
stant a in a cubic volume V = L d = Na d with peri- 
odic boundary conditions. This Hamiltonian will be de- 
noted by H . The Jacobian of the linear transformation 
(fj — > (f>(k) is \dtpj/d(f(k)\ — {aL)~ dN / 2 . The dimension- 
less partition function is 



vG 



N 

n 

3=1 

n 



d(pj 



exp(-iP 



aL d / 2 



d0(k) 



n 



2tt 



exp(-£T G ) 

1/2 



k v [r + ^(k)]5 2 / 
For the transformed system one obtains 
1 



(B.l) 



Z 



/G 



n 

. k' 



v n/d i>d/2 
2tt 



#'(k') 



exp(-iJ' G ) 
1/2 



[r + <SA'(k')K 2/d 



(B.2) 



with 1/ = (det A)" 1 / 2 ^. Eqs. (fB~T|) and pO]) define the 
integration measure f d</?(k) and f d<p'(k.'). The Gaus- 
sian free energy densities divided by fc^T are 



f' G = 



ln(27r) 
2a d 

ln(2^) 
2v' 



^^ln{[r + 5A(k)]fi 2 }(B.3) 



2L d 



^^ln{[r + ^'(k')] V ' 2/d }. 



2L' d 



(B.4) 

The correctness of the additive constant of f G can be 
checked by performing the integrations of Z G in real 
space for Kij — 0, <5A"(k) = 0, 



2^ 



2?r 
r a 2 



N/2 



n / exp 

i =1 -00 

(B-5) 

This is valid also for free boundary conditions. The ad- 
ditive con stant of f G was not correct in previous work 
[l02l 103], In order to calculate (|B.3[) we consider 

A(r , L, K id ,a) = m iN + <^A(k)]a 2 } 



- / ln{[r + ,5A(k)]a 2 } (B, 
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where the sum ^ k and the integral J k have finite cutoffs 
±7r/5 for each k a [see Eq. (|4.3ip ]. Using the integral 
representation 



lnw = dyy 1 [exp (-y) - exp (-wy)] (B.7) 



and interchanging the integration J dy with ^ k and J k 
we obtain, because of L~ d ^ k 1 = f, 1, 



Alro.I.ify.a) = dyy-'e-^y / exp{-5A(k)5 2 y} 



(B. 



L- d ^cxp{- ( 5if(k)a 2 2/ } 



Since (5A(k) is a periodic function of each component fc a 
of k the sum in (|B.8[) satisfies the Poisson identity [13, HH 

L- rf ^cxp{-<5^(k) ( i 2 2/ } 

k 



where n = (rii, r«2, ?v) and k • n = Y] t k a n a . The 
sum Y2 n runs over all integers n a , a — 1, 2, in the 
range — oo < n a < oo. This leads to the exact represen- 
tation 



A(r ,L,K lj7 a) = - J dyy-'e-^y 
o 

x / exp{-5#(k)a 2 t/ + ik • nL}. (B.10) 

Note that here the integral f. still has finite lattice cutoffs 
±7r/o. We shall evaluate A(ro, L, Kij,a) for large L/d^S> 

1 and distinguish two regimes: (i) Lr^/ 2 < 0(1), r c/ 2 o ^ 
1, and (ii) Lr^ 2 3> 1 for arbitrary fixed r^ 2 a > 0. 

In the regime (i), the large - k dependence of SK(k) 
does not matter. Therefore we may replace 5K(k) by its 
small - k form k ■ Ak and let the integration limits of 
J, go to oo. Furthermore it is useful to substitute the 
integration variable z = 47r 2 a 2 j/ / L 2 . Then we obtain 



A(r ,L,K id ,a) -» A(r , L; A) = - / cto^e-^ 2 ^ 4 ^ 



The Gaussian integral over k yields 



oo 

J exp[-k • AkL 2 z/(47r 2 ) + ik ■ nL] 



(detA)- 1 / 2 ( 4- N ) d/2 exp(-n- A^n^/z) .(B.12) 



Eqs. (jB~TT]) and (|B~T2|) lead to 

A(r ,L;A) - i- d g (r L' 2 , A) 



(B.13) 



where Go and Kd(y,A) are given by (|6.34|) and (I4.46[) . 
Thus, in the regime (i), we derive from Eqs. (|B.3|) . (|B.6I) . 
dEU, and (|BTT5|) 



f° = fb+ ^L- d g a (r L ,2 ;A) , (B.14) 



where the bulk part f G is obtained from (IB.3jl by the 
replacement L _d ^ k — > Jk- We note tnat within the 
anisotropic Gaussian model there exists no unique corre- 
lation length. This exists only for the transformed sys- 
tem with the (asymptotically isotropic) Hamiltonian H' G 
[Eqs. (|2~^0j) and (j^g]) for u' = 0, h' = 0] for which the 
corresponding result reads 



f iG _ f iG 

J — Jb '2' 



L'- d g (r L u ;A) , 



(B.15) 



with the bulk part f' G given in (|4.29|) for u' Q = 0. Now 
the parameter rg is related to the second-moment bulk 
correlation length of H lG [see (|3.4[) ] 



e; G = r - 1/2 = ^.r 1 ' 2 p' G = - 1/2 



£'1-1 

?0H 



(B.16) 



This leads to the identification of the scaling function of 
the Gaussian excess free energy density in the regime (i) 



T G > EX (£;A) = (l/2)g Q (i;A) 



(B.17) 



with i = tiL'/^) 1 ^, v = 1/2. 

In the regime (ii), A(ro, L, K^j, a) will attain an ex- 
ponential L dependence and the complete k depen- 
dence of SK(k) does matter. For simplicity we consider 



only nearest-neighbor couplings A^. 



K on a sc lat- 



tice, SK(k) = 4A"^ Q , =1 [1 — cos(afc Q )], with the second- 
moment bulk correlation length of the Gaussian model 



, 2K 

£ = a — 



ro 



1/2 



2K 

«o 



1/2 



(B.18) 

At the level of the Gaussian model there exist no Wegner 
corrections to (|B. 18|) . Using (|B.7[) we obtain, similar to 
App. A of [77[, the exact representation for (|B .6[) 



( cxp[-k • AkL 2 z/(47r 2 ) + ik ■ nL] . (B.ll) 

1*0 £ 



- 1 P ~ r oV 



A(r ,L,K,d)=d dyy e 



x { [5(oo, y)] d - [S(L/a, y)] d } , (B.19) 
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S(L/a,y) = S(oo,y)+2e- 2 y £ 7 mi/a (2y) ,(B.20) 

m— 1 

with SXoo, y) = e- 2 yi (2y) and f = r /(2K) = (~a/^) 2 
where 



Jo(rL' 2 ,A) 



dy 
y 



exp 



rL'V 



{(tt/z/)^ 2 - K d (y, A) + l} - exp(-y)] . (C.2) 



I M (2y) = (2tt)- 



d6> cos(A/6>) e 



2y cos £ 



(B.21) 



are the m odifie d Bessel functions with integer M (see 
(9.6.19) of [To3]). In the limit of large L/a for fixed f > 
only the range of y ~ 0(L/a) is relevant and only 
the m = 1 term of (|B.20[) suffices to obtain the leading 
exponential behavior. Thus we substitute 2y = zL/a in 
(|Bl9|> and use the asy mpt otic formulae for larg e L /a (see 
(9.7.7) and (9.7.1) of [To| and App. A of pjf ) 



h/d(zL/a) ~ {2-nLq/a) 



a 



hi 



I (zL/a) ~ e zi / 5 (27rzL/5)- 1 / 2 



(B.23) 



where g = (1 



r2U/2 



The maximum of the resulting 



exponential part of the integrand of (|B.19[) is at 



where z 



[fo(l + f /4)] 



-1/2 



Expanding around z 



and performing the integration over z yields for large L/a 
at arbitrary fixed f o > 



with the exponential correlation length 



e G = 

Se 



rsinh 



2^ 



(B.25) 



No universal finite-size scaling function of the Gaussian 
model can be defined in the region 3> 1 because of 

the explicit a dependence of (|B .24|) and (|B.25[) . 

Within a RG treatment of the ip 4 lattice model the 
Gaussian results can be considered as the bare one-loop 
contribution. By means of such a RG treatment at fi- 
nite lattice constant a parallel to Sect. 2 and App. A 
of the results derived above acquire the correct crit- 
ical exponents of the n = 1 universality class including 
corrections to scaling. This leads to the one-loop results 
(|10.14p - (|10.16p which are valid for arbitrary 5/£± > 0. 

For field theory with SK(k) — k 2 and a sharp cutoff 
— A < k a < A the Gaussian excess free energy density 
is given by (1/2) Aa where Aa is given by (|B.6D with a 
replaced by A -1 . From (A. 31) of App. A of [89] and a 
RG treatment at finite A we obtain (|10.27|) - (|10.31|) . 
Appendix C : Sums over higher modes 



The v' dependent finite-size part in (|C.1[) comes from the 
absence of the k' = term and is exactly cancelled by the 
corresponding logarithmic term in (|4.26p . For d > the 
function Jo (rL 12 , A) is finite for all < rL' 2 < oo and 
diverges for large rL' 2 as Jo(rL' 2 , A) ~ — \n[rL' 2 /(4ir 2 )]. 
By means of differentiation with respect to r we obtain 



(B.22) 



S m (r)=L'- d J2[r + SK'(k' 



IClJW—m 



r + 5K'(k') 



(27T) 



2m 



I m {rL ,A), (C.3) 



I m (rL' 2 ,A) = / d 2 /j/ m - 1 exp[-rL' 2 2/ /(4^ 2 )] 



x{K d (y,A)-(n/y)V 2 -l} . (C.4) 



For 2 < d < 4 the behavior of these functions for r — > 
is h(r,A) -> Ii(0,A) = finite, r 2 / 2 (r, A) - 0(r d / 2 ), 



Mr,A) ^h(r,A) - ^I 2 (r,A) 
= J (0,A) + 0{r d ' 2 ) . 



(C.5) 



For d > the behavior for r — > oo is rl±(r, A) — > —Air 2 , 
r 2 I 2 {r, A) -> -16tt 4 . For 2 < d < 4 and r > the bulk 
integral in (jC.lj) is 



ln{[r + 5K'(k')}v' 2/d } = J \n{[5K' {k')]v' I/d } 

k' k' 

+ r [[6K\k')]- 1 -2A d r d / 2 /(de),(G.6) 



Using (|B.6|) and the integral representation (|B.7[) with 
to = 7\L /2 /(4tt 2 ) we obtain from (|BTT3|) and ([6~M)l 

* £ln{[r + 5X'(k')]t/ 2 / d } = \ [ln{[r + 6K'(k')}v' 2 /*? 

k'^O 1,1 



apart from nonasymptotic corrections that depend on 
jj/2/<2 anc j van j sn f or r _^ o at fixed finite w'. The bulk 
integrals in (|C.3|) follow by differentiation with respect to 
r. 



2L' d 



In 



t ,/2/d4 7r 2 



+ 2l^ J ° (rL ' 2 ' A) ' (C - 1} 
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